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Abstract. The pseudogap Anderson impurity model provides a paradigm 
for understanding local quantum phase transitions, in this case between 
generalised fermi liquid and degenerate local moment phases. Here we develop 
a non-perturbative local moment approach to the generic asymmetric model, 
encompassing all energy scales and interaction strengths and leading thereby to 
a rich description of the problem. We investigate in particular underlying phase 
boundaries, the critical behaviour of relevant low-energy scales, and single-particle 
dynamics embodied in the local spectrum D{u>). Particular attention is given to 
the resultant universal scaling behaviour of dynamics close to the transition in 
both the GFL and LM phases, the scale-free physics characteristic of the quantum 
critical point itself, and the relation between the two. 
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1. Introduction 

The essential physics of a single magnetic impurity coupled to a fermionic host is 
embodied at its simplest in the Anderson impurity model (AIM) [1]: a correlated, 
non-degenerate impurity with local interaction U, hybridizing to a non-interacting 
host band with density of states /9host('^) (for a comprehensive review see [2]). Yet 
its simplicity is nominal — even for the conventional case of a metallic host, the 
basic model for imdcrstanding magnetic impurities in metals [1,2] and highly topical 
again in the context of quantum dots for example [3,4], or surface atoms probed by 
scanning tunneling microscopy (STM) [5]. Here the essential strong coupling (large- 
U) behaviour is that of the Kondo effect, characterised by a single low-energy Kondo 
scale ujk, manifest vividly in the Kondo resonance in the local single-particle spectrum 
D{w), and rendered no less subtle for being by now quite well understood. 

There is however one sense in which the metallic AIM is 'simple': it lacks a 
non-trivial quantum phase transition. The impurity spin is quenched ubiquitously by 
coupling to the low-energy degrees of freedom of the metallic host, the system is a 
Fermi liquid for all interaction strengths U and the Kondo scale, while exponentially 
small in strong coupling, never quite vanishes. 

That situation changes markedly if the host contains instead a pseudogap in 
the vicinity of the Fermi level [6] (cu = 0), namely /Ohost(^) oc \lu\^ with r > 
(and r = corresponding to the normal metallic AIM). This pseudogap Anderson 
impurity model (PAIM) has been studied extensively in recent years [6-20], particularly 
for large-C/ where the low-energy subspace maps onto that of the associated Kondo 
model, and from which work much progress and insight has accrued. The resultant 
depletion of host states around the Fermi level leads in general to the occurrence [6- 
20] of a quantum phase transition at a finite, r-dcpendent Uc- Above the quantum 
critical point (QCP) a degenerate local moment (LM) phase arises, with a residual, 
unquenched impurity spin. For U < C/c by contrast, and perturbatively connected to 
the non-interacting limit (U = 0), resides a 'strong coupling' or generalised Fermi 
liquid (GFL) phase in which the impurity spin is locally quenched and a Kondo 
effect manifest, reflected in a low-energy Kondo scale ojk- In the vicinity of the 
transition each phase is in fact characterised by a low-energy scale [17-20], here denoted 
generically by lu^; corresponding to the Kondo scale lu^ = luk for the GFL phase and 
its counterpart = wl [20] for the LM phase. Most importantly, and symptomatic 
of the underlying transition, the a;*'s vanish as the QCP is approached. 

Although the PAIM serves in general as a paradigm for understanding local 
(boundary) quantum phase transitions [21], many examples of specific systems with 
a power-law density of states have also been discussed. These include various zero- 
gap semiconductors [22], one-dimensional interacting systems [23], and the important 
problem of impurity moments in d-wave superconductors [17,18], notably Zn- and Li- 
doped cuprate materials. Likewise it has recently been shown [19] that the pseudogap 
Kondo model exhibits critical local moment fluctuations, and an associated destruction 
of the Kondo effect, very similar to those present at the local QCP of the Kondo lattice 
[24] that has been argued to explain the behaviour of heavy fermion metals such as 

CeCu6-a;AUa;. 

While great progress has been made in understanding the PAIM [6-20], much 
nevertheless remains to be understood; particularly regarding the quantum critical 
behaviour of dynamical properties close to and on either side of the transition, and 
precisely at the QCP itself. Since the low-energy scale is arbitrarily small (but non- 
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zero) sufficiently close to the transition, the frequency {u>) dependence of dynamics 
should scale in terms of it — i.e. exhibit miivcrsal scaling as a function of uj/lu^ — 
with resultant scaling behaviour that will depend on the particular phase considered, 
GFL or LM. What precisely is that behaviour for the selected dynamical property? 
And how from it docs (or can) one imderstand the scale-free physics characteristic of 
the QCP itself, where = identically? These are tough issues, especially from an 
analytical perspective, and some key aspects of which have recently been uncovered 
via numerical rcnormalization group (NRG) calculations [16,18,19], the local moment 
approach [15,16,20] and, for the multichannel pseudogap Kondo model, a dynamic 
large- A'' approach [17]. Single-particle spectra embodied in the impurity D{u}) provide 
an obvious example of dynamics [12, 14-20]; and an important one given e.g. the direct 
manifestation in the GFL phase of low-energy Kondo physics embodied in the Kondo 
resonance, as well as their more general relevance to STM and related experimental 
techniques. 

In this paper we develop a local moment approach (LMA) [25-30,15,20] to the 
PAIM at T = 0, for the general case of the particle- hole asymmetric model (noting that 
the corresponding Kondo model thus contains both exchange and potential scattering) . 
Single-particle dynamics form a central focus of the theory albeit, importantly, that 
its intrinsic formulation also enables phase boundaries between the GFL and LM 
phases to be determined, as well as the critical behaviour of associated low-energy 
scales [15,20]. Although developed originally to describe metallic Anderson impurity 
models, at T = [25-27], finite temperatures [28] and magnetic fields [29], the approach 
is quite general; readily extended for example to encompass lattice-based models 
within the framework of dynamical mean-field theory, such as the periodic Anderson 
lattice model [30]. It has also been applied previously to and has led to a number of 
new predictions for, the symmetric limit of the PAIM [15,20]; rather successfully, as 
judged by subsequent comparison to NRG calculations [16]. The asymmetric PAIM 
considered here is however another matter: the symmetric model, recovered though 
it is as a particular limit of the present work, has long been argued [13,17,18] to be 
a somewhat special case, in some sense divorced from generic asymmetric behaviour 
(although our own caveats on the matter will become clear in due course). 

The LMA, while certainly approximate, is also well suited to the present 
problem. Its intrinsically non-perturbative nature enables it to handle strong electron 
correlations and the attendant, central low-energy scales; while at the same time 
permitting access to dynamics on all energy scales, whether it be on the lu^ scales 
relevant to universal scaling behaviour, or scales on the order of a; ~ i7 relevant to 
the high-energy Hubbard satellites. Importantly too, and in contrast to difficulties 
inherent to a number of other potential theoretical approaches, the method handles 
both the GFL and LM phases on an essentially equal footing (the LM phase, for 
example, is certainly not a trivially decoupled/isolated spin). 

The paper is organised as follows. After the requisite background, we consider in 
§2.1 an obvious question: since the GFL phase is pertubatively connected to the non- 
interacting limit, satisfying as such the Luttinger theorem [31], what is the analogue of 
the resultant Friedel sum rule? That is entirely familiar for the conventional metallic 
AIM (r = 0) [2,32], and its counterpart for the PAIM plays a central role in the 
present work. In §3 another important issue is considered: what can be deduced on 
general, minimalist grounds about the behaviour of the scaling spectra close to the 
quantum phase transition; in both the GFL and LM phases and at the QCP itself? 
Such considerations are revealing, and dictate constraints to be met by any credible 
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approximate theory. They show moreoever that a variety of low-energy scales in 
addition to a;* enter the scaling behaviour of dynamics the quasiparticlc weight Z 
and renormalised level ej* for the GFL phase, together with the latter's spin-dependent 
counterparts for the LM phase. Each vanishes as the transition is approached, with 
distinct but related critical exponents; but nonetheless enter dynamics in such a way 
that universal scaling remains determined solely by the key low-energy scale w*. 

The LMA itself is developed in §4, focusing on the notion of symmetry restoration 
that is central to the approach [15,20,25-30], together with what likewise proves 
important in dealing with the asymmetric PAIM, namely self-consistent satisfaction 
of the Priedel sum rule/Luttinger theorem. Results arising from the approach are 
presented in §s 5-7, beginning with the critical behaviour of relevant low-energy scales 
(§5.1) and a determination of representative phase diagrams for the problem (§5.2). 
Single-particle dynamics are considered in §6, on all energy scales whether universal or 
non-universal; but with a dominant focus on scaling spectra in the GFL and LM phases 
— on all lv/lo^ scales — as well as the scale- free dynamics arising at the QCP. Finally, 
we consider in §7 the regime of small r <C 1, obtaining explicit analytic results for the 
critical behaviour of low-energy scales, resultant phase boundaries and single-particle 
scaling spectra. The paper concludes with a brief summary /outlook. 

2. Background 

With the Fermi level taken as the energy origin, the Hamiltonian for the spin-| AIM 
[1,2] is given by 



The first term describes the non-interacting host with density of states Phost('^) = 
J2k ~ ^k), while the second refers to the correlated impurity with bare site-energy 
Ci and on-site interaction U. The third term is the one-electron host-impurity coupling 
embodied in Vifc, taken for convenience as Vn, = V (although that is not required). 
The particle-hole asymmetry of the model is specified by the asymmetry parameter 



the impurity charge n = '^^{fiia) = 1 for all U. In strong coupling, where the relevant 
low-energy sector of the AIM maps under a Schrieffer- Wolff transformation onto a 
Kondo model containing both exchange (J) and potential (K) scattering contributions 
[2], the asymmetry rj is equivalently the ratio of the potential/exchange scattering 
strengths: r] = 1 - 2\ei\/U = K/J (see e.g. [27]). 

We focus here on the impurity Green function G(a;)(<-> G{t) = — i(T (cio-(i)cj^))), 
and hence the local spectrum D{lj) — -'^sgn{uj)lm.G{(jj). In addition to intrinsic 
interest in single-particle dynamics per se, this also provides a direct route to both 
the location of phase boundaries between the GFL and LM states [15] and the critical 
behaviour of associated low-energy scales, as will be shown in §s 4,5. 

G{lo) is conventionally expressed as 




(2.1) 




■i = --f and 



(2.2) 



G{oj) = [lu+ -€i- A{uj) - E(w)] 



-1 



(2.3) 
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with u}~^ = u) + iO+sgn(a)) and S(a') = — isgn{u)T} (uj) the usual single 

sclf-cncrgy; the corresponding non-interacting Green function is denoted by gi^uj). 
All effects of host-impurity coupling at one-electron level are embodied in the 
hybridization function A{uj) = Y^^^iki^'^ ~ ^k]~^ (= Ar — isgn(a;)Ai(w)). For 
the PAIM considered, Ai{uj) = 7rF^/9host('^) is given explicitly by Ai{uj) — 
Ao{\uj\/Aoy0{D — with Aq the hybridization strength and D the bandwidth 
{0{x) is the unit step function). That this form is naturally simplified is of course 
immaterial to the relevant low-energy physics of the problem in the scaling regime of 
primary interest. We take Aq = 1 as the basic energy unit in this paper, i.e. 

Ai{lu) :^\uj\''e{D -\uj\) (2.4) 

with r > for the pseudogap AIM (and the limit r = corresponding to the normal 
metallic model). The real part Aji{u) follows from Hilbert transformation and, for 
later reference, has the low-o; behaviour (for all r > 0) 



Anioj) = -sgn(a;) J f3ir)\u;\^ + + O 

I 7r(r ~ 1) D 



D 



(2.5) 



with (3{r) = tan(^r) used throughout. Notice from this that the wide-band limit 
Z) ^ 00, so frequently employed for the metallic AIM (see eg [2]), can be taken for 
r < 1, whereupon AR(a;) = —sgn{uj)(3{r)\uj\'' for all lo. We consider < r < 1 in this 
paper, which range is assumed unless stated to the contrary; and for that reason will 
employ the wide-band limit. This is simply a matter of convenience, is again irrelevant 
to the key low-energy behaviour of the system and can be relaxed trivially. 

The impurity charge n is of course related to the local spectrum by n = 
2 J^^dw D{u)). The 'excess' charge nimp will also prove important to the present 
work, for reasons given in §2.1. This is defined as the change in the number of electrons 
of the entire system due to addition of the impurity to the host, and is given generally 

by 



J — <x 



du Ap{uj). (2.6) 



Here Ap{uj) (= p(a;) — phost{<^)) is the excess density of states, itself given by (see e.g. 
[2]) 



Ap{lj) = — — sgn(a;)Im \ G{lo) 
n 



dA{w) 



dui 



(2.7) 



such that riimp is related to n by 



ni^p = n- h^n J du G{u)^^^. (2.8) 

From the above, for specified gap index r, the PAIM is characterised by any two 
of the 'bare' parameters U, Cj and r] = 1 + 2ei/U (together if appropriate with the 
bandwidth D). In practice, as now explained, we choose to consider ry and one or other 
of U and ei as the basic parameters. For the r = metallic AIM we have recently 
considered [26,27] the scaling behaviour of -D(u;) arising in the strong coupling Kondo 
regime, where Ci < such that |ei| ^ 1 and U — \e\\ ^ 1. Here the Kondo scale 
is exponentially small. For fixed asymmetry 77 [27], D{u}) then exhibits scaling in 
terms of a) = w/wk- That is, on fixing rj = 1 — 2\ti\/U and progressively increasing 
the interaction U, D{u) is universally dependent on a) = w/wk alone, and otherwise 
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independent of the bare model parameters (which themselves determine the Kondo 
scale). The important point here is that this arises for fixed asymmetry t/ = K/J 
(and in this case only for fixed 77, i.e. universal scaling does not arise upon increasing 
U for fixed |ei| » 1 [27]). In consequence, a family of scaling spectra arises, one for 
each asymmetry 77. This is not widely appreciated (although in the related context of 
d.c. transport it was understood early by Wilson [33]), but that it must arise is nearly 
obvious. For if the scaling spectra for the generic asymmetric AIM were 7?-independent, 
they would coincide with that for the p-h symmetric model 77 = (= K/J) and as 
such be strictly symmetric in u> about the Fermi level w = 0; the implausibility of 
which is at least intuitive. 

For the pseudogap model considered here, the quantum phase transition between 
GFL and LM phases is naturally characterised by a low-energy scale — the Kondo 
scale wk for the GFL phase and its counterpart wl [20] (discussed in §6ff) for the 
LM phase that vanishes as U Uc{r)± and the transition is approached. In 
consequence the spectra likewise exhibit scaling in terms of d; = w/w* in the vicinity 
of the transition. As discussed in §6ff, we again find this scaling to be ry-dependent. 
For this reason, and to enable continuity as r ^ to the above scaling behaviour 
for the metallic AIM, we choose to regard ry and one or other of U and Ci as the 
basic model parameters. The minimal range of 77 that need be considered is also 
readily established, for under a general p-h transformation of the Hamiltonian it is 
straightforward to show that ei — [ei + U] and -D(a;; Ci) = D{—lo; — [ei -|- U]), hence 
n(ei) = 2 — n(— [ci + [/]) for the impurity charge (and likewise for the excess charge 
n\aii>)- From this it follows that only ei > —\U need be considered, corresponding to 
77 > and which range is assumed from now on. 

To capture the underlying phase transition intrinsic to the PAIM, an inherently 
non-perturbativc approach is obviously required. That poses severe and, to 
our knowledge, insurmountable problems for any approach based directly on the 
conventional single self-energy Y,{uj). But such an approach is not obligatory, S(w) 
merely being defined by the Dyson equation implicit in equation (2.3) for G{uj). For 
the doubly degenerate LM phase in fact, a two-self-energy (TSE) description would 
be physically far more natural (as is obvious by considering e.g. the atomic limit, an 
extreme case of the LM phase). Here the rotationally invariant G{w) is expressed as 



in terms of self-energies Scr(w) (= S^(a;) — isgn(a;)Sj,(a;)) that are in general distinct 
and spin-dependent. The TSE description is in general a necessity and not a luxury 
for the LM phase. It is precisely this framework that underlies the LMA [25-27] where, 
as explained further in §4, it is employed for both the LM and GFL phases and in 
consequence is able to capture simultaneously both phases and hence the transition 
between them. The conventional single self-energy can of course be obtained as a 
byproduct of the TSE description. It follows from direct comparison of equations 
(2.3,2.9,2.10) and is given for later use by 



G{u;) = \[G^{uj) + G^{uj)] 



(2.9) 



with the G„{uj) given by 



G,(w) = [u;+ - e; - Aiuj) - E,(a;)] 



-1 



(2.10) 



S(a;) 




i(|]T(^)-Ex(c.)) 



(2.11) 



with g{uj) = [uj 
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Before turning to the LMA itself we consider two important questions. What 

can be deduced on general grounds about the behaviour of the GFL phase? That is 
considered in §2.1. Second, as considered in §3, what may be inferred via general 
scaling arguments about the behaviour of the scaling spectra close to the phase 
transition; in both the GFL and LM phases and at the quantum critical point (QCP) 
itself, where the low-energy scale w* = identically? 



2.1. GFL phase: Friedel sum rule 

The GFL phase is perturbatively connected to the non-interacting limit, which in 
essence defines that phase. This adiabatic continuity requires that the Luttinger 
integral vanish [2,31], i.e. that 

doj G{u) = 0. (2.12) 



= lm f 

J —cx 



From this, using only S'(0) = 0, the Friedel sum rule follows directly. Given explicitly 

by 

nimp = 1 - -tan-i ( -r^ ) (2.13) 



it relates the excess charge riimp to the 'renormalised level' e*, itself defined by 

er = ei + SR(0). (2.14) 

For the r = metallic AIM, where Ai(0) = Aq (=1), this gives the Friedel sum rule 
in the famihar form e* = tan[|(l - nj^p)] [2]. For the PAIM by contrast, Ai(0) = 0; 
and in consequence nimp is restricted to the integer values 

'o - t^n (2-15) 
2 : ef < ^ ' 

according to whether the renormalised level lies above or below the Fermi level. 
This behaviour, which we reiterate is a consequence of adiabatic continuity and the 
Luttinger theorem, is physically natural: one expects intuitively that the change in 
number of electrons due to addition of the impurity should be integral. It is in fact 
the more familiar metallic AIM, where from equation (2.13) above riimp is in general 
non-integral, that is the exception; reflecting in that case the existence of a finite 
density of conduction band states arbitrarily close to the Fermi level. 

The above conclusions may be checked order-by-order in perturbation theory (FT) 
about the non-interacting limit (as we have confirmed explicitly at simple first and 
second order levels) : the resultant ni„ip calculated from spectral integration (equation 
(2.6)) indeed satisfies equation (2.15), and Jl = 0. However the Friedel sum rule and 
hence Luttinger theorem is not satsified by approximate theories in general, save for 
the p-h symmetric case where it is guaranteed by symmetry (and for which Ui^p = 1 
regardless of whether the GFL or LM phase is considered). For the metallic AIM 
that is well known to be the case for the non-crossing approximation (NCA) [34] , and 
is true even for second order PT employing Hartree renormalised propagators [35], 
usually but loosely referred to as second order PT: the resultant nimp determined by 
spectral integration does not in general concur with that implied by the Friedel sum 
rule (whether for the normal metallic model or the PAIM). It is however obviously 
desirable that the Friedel sum rule be satisfied if possible, and in §4.1 we incorporate 
it straightforwardly within the LMA. 
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3. Scales and scaling: general considerations 

Close to and on either side of the quantum phase transition intrinsic to the PAIM, the 
problem is characterised by a single low-energy scale denoted generically by u* that 
vanishes as the transition is approached, U — > Uc{r)± [17-20]. Given this we now use 
general scaling arguments to make a number of deductions about the scaling spectra, 
in both the GFL and LM phases and at the quantum critical point (QCP) itself 
(where oj^ = 0). We take it for granted in the following that the QCP is continuously 
accessible from either the GFL or LM phases as ^0. Although natural, implicitly 
assumed in NRG studies [18,19] and found also to arise within the LMA [20], this 
remains strictly an assumption (without which essentially no deductions can be made 
about the QCP). 

The spectrum is denoted in this section by D{U;u)) with the {/-dependence 

temporarily explicit (and the rj = 1 + 2ei//7 dependence implicit). As the transition 
is approached, u = \1 — U/Uc{r)\ and the low-energy scale vanishes, as 

= u" (3.1) 

with exponent a. D{U;u}) can now be expressed generally in the scaling form 
ttD{U; u) = m^"''5'q(cj/m") in terms of two exponents a and h (and with a = GFL or 
LM denoting the phase), i.e. as 

7Tu;tD{U;uj) = -^jL.'/iu,) (3.2) 

with the exponent a eliminated and the w-dependence encoded solely in u = w/w*. 
Equation (3.2) simply embodies the universal scaling behaviour of the single-particle 
spectrum close to the transition. Since the QCP must be 'scale-free' (independent 
of w*), it follows from equation (3.2) that the leading x ^ oo behaviour of ^'^(a;) is 
independent of a and given by 

^a{Lo/u.) "/"Ar°° C(r, r;)(|w|M)-'' (3.3) 

with C(r, rj) a constant; such that precisely at the QCP (w* = identically) 

nD{U,;w) = C{r,r])\u>\-K (3.4) 

The exponent b thus governs the w-dependence of the QCP spectrum, and is 
independent of the phase a. In physical terms equation (3.3) embodies the fact that 
the large w asymptotic 'tail' behaviour of the scaling spectra in either phase coincides 
with that of the QCP. More generally of course, the (D-dependence of the scaling 
spectra ^'a(tj) will differ for each phase. 

For the p-h symmetric PAIM (r/ = 0) it is readily seen that the exponent b = r, 
using equation (3.2) and the fact [14] that the leading low-w behaviour of D{U; co) in 
the GFL phase is entirely unrenormalised from the non-interacting limit and given by 

7rD(f/;a;) ''^i^°cos2(|r)|a;|-'' : r] = (3.5) 

We add that the resultant QCP behaviour TrD{Uc;u)) = C(r; 0)|a;|~'' is precisely as 
found in an LMA study of the symmetric PAIM [20], and that the exponent 6 = r is 
likewise consistent with NRG calculations for the GFL phase [16]. 

The behaviour equation (3.5) is exact [14-16], arising because S(a;) vanishes more 
rapidly than the hybridization (oc [^1'"), which itself reflects adiabatic continuity to the 
non-interacting limit. That is manifest also in the limiting low-frequency 'quasiparticle 
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form' for the impurity Green function in the GFL phase. This follows from a low-w 
expansion of the single self-energy 

~ SR(0) - (^1 - 1^ a; (3.6) 

with Z = [1 — {d'£^{u))/dLv)u;=o]~^ the quasiparticle weight which likewise vanishes as 
U Uc{r)—; and with quasiparticle damping embodied in E'(a;) neglected as usual 
[2], as is readily justified perturbativcly. For the symmetric PAIM, where e; = —U/2 
and the impurity charge n = 1, I]^(0) is given trivially by the Hartree contribution 
^Un = \U, whence the renormalised level e* = ei + S^(0) vanishes by symmetry. 
From equations (2.3) and (3.6) the resultant quasiparticle behaviour is then given in 
scaling form by 

7ra;:D(C/;w) ~ 7— , Jj^,-, ^n — ■ V = (3-7) 

* ^ ^ {kiuj + sgn{tj)(3{r)\u)\'-)^ + lu}]^"- ' ^ ' 

from which equation (3.5) is recovered as |d)| — > 0; and where the Kondo scale and 
quasiparticle weight Z are related by 

= [fciZ]T^ (3.8) 

(which relation we add is indeed found within the LMA [15]). 

For the general asymmetric PAIM {r\ ^ 0) we know of no simple argument to 
infer the exponent 6, but NRG calculations arc again consistent with b = r [18] (as 
one might expect from continuity to the symmetric limit). For the asymmetric model 
the renormalised level e* will in general be non-vanishing, and the quasiparticle form 
appropriate to the GFL phase is then given by 

la)!'" 

(fciw - {e*/uji) + iigii(LU)p(r)\Lu\'^y + 

Granted 6 = r it follows from this that the renormalised level e?" must also vanish as 
U Uc{r)— and must do so at least as rapidly as ul, i.e. 

|en~w"' (3.10) 

with exponent cii > ar. The GFL scaling spectrum ^'GFL(a)) = t^^ID{U:ll!) is of 
course obtained by considering in general finite w = w/w*, in the formal limit — > 0. 
From equation (3.9) there are then two a priori possibilities for its \ow-Q behaviour, 
(a) Either ai > ar and thus f* /col 0; whence from equation (3.9) the leading low-d; 
behaviour of the asymmetric GFL scaling spectrum would coincide with that of the 
symmetric limit equation (3.7), namely loID{U;ijj) oc \oj\~^ . Or (b) ai = ar and the 
renormalised level vanishes precisely as col, 

u;. = [\e*\/k]i (3.11) 

with k a constant. In this case the leading low-w behaviour of the scaling spectrum 
is ^'gfl(w) C)c Id'l'', and the quasiparticle form equation (3.9) implies an asymmetric 
double-peak spectrum straddling the Fermi level. It is precisely this characteristic 
behaviour that is observed in NRG calculations [18] (see also §6.1 and figure 11 below), 
from which we conclude that and the renormalised level should indeed be related 
by equation (3.11). NRG calculations for | < r < 1 [18] are moreover consistent with 
an exponent of a = i for the vanishing of the low-energy scale uj^ (oc T* in the notation 
of [18]); implying in turn that the exponent ai = 1. The behaviour just inferred will 
be considered within the LMA in §5,6. We would however add that it gives by no 
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means the whole story, for example the quasiparticle form is by definition confined 
to asymptotically low frequencies Cj and essentially nothing can be deduced from it 
regarding the general form of \1/gfl('^). 

The latter remarks have focused on the GFL phase. But the renormalised level e* 
will in general be non-zero in the LM phase as well, with the leading low-uj behaviour 
of the spectrum thus given from equation (2.3) by nD{U;Lu) oc |w|'"/|ej*p (assuming 
only that T,{ui) vanishes no less rapidly than the hybridization as |w| 0). For this 
behaviour to be part of the scaling spectrum in the LM regime, and because the 
exponent b (= r) of equation (3.2) is independent of the phase, it follows that the 
low-energy scale is again related to the renormalised level by equation (3.11); and 
hence that |e*| must vanish as the transition is approached from either phase. This 
too will be considered within the LMA in §5.1. 

4. Local moment approach 

In this section we develop a local moment approach to the PAIM for arbitrary 
asymmetry 77. For the particular case of p-h symmetry (rj = 0), it reduces to the 
approach given hitherto in [15]. Further details can be found e.g. in [15,27]. 

4.1. Basics and symmetry restoration 

There are three essential elements to the LMA. First that local moments ('/i'), regarded 
as the primary effect of electron interactions [1], are introduced explicitly from the 
outset. The starting point is thus simple static mean-field (MF) i.e. mirestrieted 
Hartree-Fock (HF), containing in general two, degenerate broken symmetry MF states. 
Notwithstanding the severe limitations of MF by itself (see e.g. [15,25,27]), it can 
nonetheless be used as the starting point for a successful many-body theory [25- 
30,15,20]. To that end the LMA employs the two-self-energy description mentioned 
in §2, that follows naturally from the underlying two MF saddle points. The local 
propagator is thus expressed in the form equations (2.9,10), with the two self- 
energies Sct(w) = S(j[{^o-}] (f =T, i) built diagrammatically from the MF propagators 
(themselves denoted by Qr^{ui)). The self-energies are conveniently separated as 

f:,{uj) = t1 + j:^{uj) (4.1) 

into a purely static contribution (E^, figure 1 below) that alone would be retained 
at pure MF level, together with the crucial dynamical piece T,„{ui) that contains in 
particular the correlated electron spin-flip physics intrinsic to both the GFL and LM 
phases. 

For the GFL phase the third key notion behind the LMA, as now discussed 
for the PAIM, is that of symmetry restoration [26-30]: self-consistent restoration of 
the symmetry broken at pure MF level, and hence recovery in particular of the low- 
frequency quasiparticle form (§3) for G{u)) that is symptomatic of adiabatic continuity 
to the non- interacting limit. This is given explicitly by 

(4.2) 

(from which equation (3.9) above follows directly). Equation (4.2) reflects the fact 

that as \uj\ the single self-energy S^(cj) decays to zero more rapidly than the 
hybridization (oc |a;|''), so that neglect of quasiparticle damping is justified and the 
leading low-w behaviour of G{oj) amounts in general to a simple renormalisation of 



— -e* + sgn{uj)l3{r)\uj\'' + isgn(w)|a;|' 
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the non- interacting limit (§3). The question then is under what conditions on the 

two-self-energies {'^^■{uj)} does such behaviour arise? In parallel with the analysis of 
[27] this is readily answered via a simple low-o; expansion of the So-('^)'s, namely 



S^(a;)~S^(0)- 



(J (4.3) 



for the real parts, with = [1 — {dT,^{u)) /du})u,=o\ ^ thus defined; together with 

f:l{uj)^El{u;)r.aM"' (4.4) 
for the imaginary part, with m > r. This form is in fact guaranteed from the diagrams 
(§4.3) for T,„{u) = T,[{Qcr}] with m = 2 + 3r (which concurs with simple perturbation 
theory in U). Equations (4.3) & (4.4) are now used in equations (2.11) & (2.3) to 

determine the resultant \ow-lu behaviour of E(cj) and hence G{lu). From this one finds 
simply that the requisite condition for the quasiparticle form equation (4.2) to arise is 

tf{uj = 0) = Ef (w = 0). (4.5) 

If equation (4.5) is not satisfied, then the resultant leading low-w behaviour of S^(a;) 
is oc \u)\^ which, vanishing as rapidly as the hybridization, vitiates the quasiparticle 
form equation (4.2). 

Equation (4.5) is of course the familiar symmetry restoration condition for 
the GFL phase likewise found for the r — metallic AIM [26,27], and for the 
symmetric PAIM [15] (where it reduces to S^(0) = independently of spin cr, due 
to p-h symmetry). As in previous work [25-30,15,20], it amounts in practice (see 
also §4.2) to a self-consistency condition for the local moment |/i| entering the MP 
propagators {^cr(w)} from which the two self-energies {Ect(i^)} are diagrammatically 
constructed. This generates directly [25,27] a low-energy scale, seen as a sharp 
resonance (at u) = Um.) in the spectrum of transverse spin excitations Im 11"' (to) 
(§4.2); which is non-zero throughout the GFL phase and sets the timescale t ~ h/oj,rn 
for symmetry restoration and hence local 'Kondo singlet' formation. This is the Kondo 
scale, in terms of which single-particle dynamics are found to scale in §6; and with 
Um = = ivk in the notation employed above (low-energy scales are of course in 
general defined only up to an arbitrary constant). 

Imposition of symmetry restoration as a self-consistent constraint thus again 
underlies the LMA for the GFL phase. If equation (4.5) is satisfied then (i) the 
requisite low-energy quasiparticle form is recovered, as above; (ii) all self-energies 
coincide at the Fermi level, namely S^(0) = E^(0) for either a (whence the 
renormalised level equation (2.14) is given equivalently by ef = ei 4- S^(0)); and 
(iii) the quasiparticle weight Z is related to the {Zc} by Z^^ = ^{Z^^ + Z^^). It is 
moreover the limits of solution to the symmetry restoration condition that delineate 
the phase boundaries between the GFL and LM phases (§5.2), as for the symmetric 
PAIM considered in [15] (for U > Uc{r) in the doubly degenerate LM phase itself, 
discussed explicitly in §4.3, symmetry is not of course restored). 

While solution of the symmetry restoration condition equation (4.5) ensures 
correct recovery of the quasiparticle form for the GFL phase, it does not in general lead 
to the Friedel sum rule (equation (2.15)) being satisfied. We have recently developed a 
strategy for metallic AIMs [27] under which both symmetry restoration and the Friedel 
sum rule may be simultaneously satisfied within a two-self-energy framework. That 
approach may be modified straightforwardly to encompass the PAIM, as now sketched 
(and regardless of specific details of the T,a{^^) = ^allQa}] which are considered 
separately in §4.2). 
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Figure 1. Static bubble diagram contribution to So-(t^)- The wavy line denotes 
U, and the solid line the broken symmetry MF propagator. 



4.1.1. Algorithm The two-self-energies are functionals of the MF propagators 
{Gaiuj)} [27], given by 

g^{uj) = [10+ -ei + ax^ M^'T' (4.6) 
where x = |C/|/i| and e; differ in general from their MF values Xq = ^U\fj,o\ and 
eio = 6; + |C/no (with |^ol and no the local moment and charge obtained self- 
consistently at pure MF level). The {Scr(w)} thus naturally depend on both Cj and x: 
Eo-(w) = Scr(w; Ci, x); and, via equations (2.9) & (2.10), so too does G{uj) — G{ui; ej, x) 
and hence Uimp = riimp{ei, x). With such dependence made temporarily explicit, the 
symmetry restoration condition and Priedel sum rule (§2.1) for the GFL phase read 

t,f{0;eux) = t.f{0;eux) (4.7) 

and 

: eliei,x) > 



rii, 



2 : e*{eux)<0 ' ^^'^^ 

These equations are together sufficient to determine both x and ei. In practice, for 
any given U, it is natural to solve equations (4.7) & (4.8) by fixing Ci and determining 
Si, in the following way: (i) For given a solve the symmetry restoration condition 
equation (4.7) for x = ■^U\n\ (and hence (ii) The self-energies {Scr(^^; ei, .x)} arc 
now completely determined, and the bare Ci entering G{u)) (equations (2.9) &(2.10)) is 
then varied until the Friedel sum rule equation (4.8) is satisfied. All quantities are now 
known without looping and the procedure is a rapid numerical exercise. Alternatively 
and equivalently, if one wishes to specify a bare rj = 1 + 2gi/[/ (or ei) from the outset, 
the procedure may be repeated by varying ei until equation (4.8) is satisfied. 

Further, following [27] it is straightforward to show that Ci —a corresponds 
to the particle-hole transformation mentioned in §2, under which Ci — » — [ej -I- U] (i.e. 
77 —7]), nimp — > 2 — nimp and D{uj) D{—uj). As for the metallic AIM [27] we can 
thus focus solely on e; > 0. The case of e; = is of course the particle-hole symmetric 
model [15,25] where the Luttinger integral, equation (2.12), vanishes by symmetry, 
and riimp = 1 for all U; in this case step (ii) above then drops out of the problem and 
the strategy is as presented hitherto in [15]. 

4-2. Self- energies 

We now specify the particular approximation employed for the dynamical contribution 
Sct(w) to the two-self-energies. Equation (4.1) for ^^(a;) is given explicitly by 

t,{uj) = lU{n-a\fi\)+^a{co) (4.9) 
where the first term is the static bubble diagram (figure 1), S° = 
U J_^du) D^^{u;ei,x) with \p\ = |^(ei,a;)| and n = n{ei,x) given by 

l/ij = ^cr / do; D°{u>;ei,x) (4.10a) 

„ J —00 
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Figure 2. Principal contribution to the dynamical £0^(0;). Same notation as 
in figure 1; the transverse spin polarization propagator (see text) is shown as 
hatched. 



(4.10b) 



in terms of the MF spectral density D'^{llj) = —^sgn{Lu)lmQcr{Lo). The class of 
diagrams retained in practice for T,„{co) = I^criiGa}] remains unchanged from previous 
work, and captures the key spin-flip dynamics inherent to both the GFL and LM 
phases. For a full discussion of the diagrams and their physical content the reader is 
referred to [15,25,27]. Shown in figure 2, they translate to 

S,(a;) = r ^ g_„{uj - a;i)n-<^<^(a;i). (4.11) 

Here Il°'~"'{uj) is the transverse spin polarization propagator given by an RPA-like 
particle-hole ladder sum in the transverse spin channel, 

n'"-'"(cj) = °n'^-'^(t^)[i - [/°n'"-'"(w)]-i (4.12) 

where ^Il'^~'^{ui) is the bare particle-hole bubble, itself expressed in terms of the broken 
symmetry MF propagators {Qa} and satisfying the Hilbert transform 

dwi Im°n+~(a;i)sgn(a;i) 



°n+-(a;) 



(4.13) 



We add futher that since [15,27] Tl+-{uj) = (and likewise for the °n'^-'^(a;)), 

only 11"' {ui) need be considered henceforth. 

The functional dependence of the {'Scri'^)} on the MF propagators, for both the 
GFL and LM phases, remains as detailed in [15] for the symmetric model. For the 
GFL phase in particular, 11"' (w) obeys the same Hilbert transform equation (4.13) 
as °n"'"~ and equation (4.11) may thus be written as 



where 



r ^ Imn+-(aa;i) [e{i^^)gz^{uj, + to) + 9i-LUi)gl^iuj, + uj)] (4.14) 

J —CXI 



f 

J — c 



l?O(a;i)0(±u;i) 



a; — Wi ± iO"'" 

denote the one-sided Hilbert transforms such that Qa (w) 



(4.15) 



g+{u^) + g-{uj). 
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Stability The Hilbert transform equation (4.13) for 11+ (w) implies the stabihty 
condition TrRell"' (w = 0) = dwilmll^ (a;i)/|wi| > 0. From equation (4.12), 
together with Im'^II"' {u = 0) = 0, it foUows that this positivity is guaranteed only if 
< U'Re^II"' {co = 0) < 1. An explicit expression for Re°n+~(0) is however readily 
obtained [15,27], 

C/Re°n+-(w = 0) = M^^^ (4.16) 

ImI 

with \jl{ei,x)\ given by equation (4.10a). The pure MF local moment, denoted by 
lHol = |/io(ei)|, is itself given as the self-consistent solution of = \p,{ei,x = ^U\iJ,\)\. 
Equation (4.16) thus shows that the requisite stability is preserved for [/i] > |/j,ol, and 
that the boundary to stability (when t/Rell"' (0; ei,x) = 1) is precisely the pure MF 
local moment |/uo|. 

When = |/io(ei)| it follows from equation (4.16) that Imll^ (w) (equation 
(4.12)) contains a zero-frequency pole. This is of course the essential signature of the 
doubly degenerate LM phase, considered explicitly in the following section. In the 
GFL phase by contrast, = |/io| is determined by symmetry restoration (step 

(i) of §4.1.1); with S\ii\ > correctly found in practice, as required for stability. In 
consequence Imll"' (ui) contains a sharp resonance occurring at a non-zero frequency 
(jj = aJm(= ^*), the low-energy Kondo scale for the GFL phase. This is all entirely in 
keeping with the general picture presented hitherto for the symmetric PAIM, see e.g. 
figure 4 of [15]. As U is increased for any fixed asymmetry r] in the GFL phase, Um 
progressively diminishes ((5|/i| 0) and vanishes at a critical U = Uc{r, rj), where the 
resonance becomes an isolated pole at w = ((5|/x| = 0). This signals the transition to 
the LM phase, to which we now turn. 



4.3. LM phase 

Throughout the LM phase = |/xo(ei)| and hence C/Ren+~(0) = 1. In consequence 
11"' (u) contains an a; = spin-flip pole (of weight Q), together with a continuum 
contribution ^11"' (w) [15], specifically 

n+-(c^) = + 'n+-(a;) (4.17) 

with Q = [U^{dRc"U+- {uj)/du;)^=or^ > 0. lmSn+-(cc>) is given by equation (4.12) 
with 11"' ^11"* , and the real/imaginary parts of ^11^ (w) are again related by the 
Hilbert transform equation (4.13). The basic form for Eo-('^) in the LM phase then 
follows using equation (4.17) in equation (4.11), namely 

E,(a;) = QU^gZ^iio) +^^,(0;). (4.18) 

where ^'E„{u}) is given explicitly by equation (4.14) with 11"' ^H"*"". 

For any given U and 77 (= l+2ei/U) in the LM phase U > Uc{r, rj), Ci is determined 
by requiring that the change in number of electrons due to addition of the impurity 
{nimp, equation (2.6)) is integral, as for the GFL phase discussed in §2.1. Here however 
we require that nimp = 1, i-e. 



2 r 

= -Im / du G{uj) 



1 - 



duj 



= 1 (4.19) 



in common with the local moment regime of the atomic limit (ei < and ei + U > 0) 
to which we expect the LM phase generally to be connected perturbatively in V. 
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Satisfying equation (4.19) within the approach presented thus determines Cj in a 

manner directly analogous to use of the Luttinger theorem to determine e; for the 
GFL phase as described in §4.1.1 (and in fact nimp = 1 is found in practice to be the 
only possible integer solution for nimp in the LM phase). 

Symmetry is not of course restored in the doubly degenerate LM phase, i.e. 
E^(w = 0) ^ t^{Lo = 0). The difference E^(0) - Sj^(O) is however found to be 

sign definite for all U > Uc, such that as U is decreased towards Uc, I'^^iO) — ^^{0)\ 
progressively diminishes and vanishes as the transition is approached, U Uc+- In 
this way the phase boundary between LM and GFL phases may be established, now 
coming from the LM side of the transition U > Uc- And the resultant Uc = Uc{r,r]) 
is correctly found to coincide precisely with that obtained coming from the GFL side, 
U < Uc, as the limit of solution of the symmetry restoration condition equation (4.5) 
for the GFL phase (as indeed is to be expected since as J7 — » Uc+, the resultant 
vanishing of |E|'(0) — E|'(0)| corresponds to the 'onset' of symmetry restoration). 

5. Scales, exponents and phase boundaries 

We turn now to results arising from the LMA specified above. Single-particle dynamics 
will be analysed explicitly in §6fF. Here the behaviour of the relevant low-energy scales 
in the vicinity of the quantum critical point is first considered, followed in §5.2 by 
resultant phase diagrams in the {U,r]) (or equivalently {U, ei)) plane. 



Our aim here is simply to establish the critical behaviour of the relevant low- 
energy scales, specifically cj* (in terms of which single-particle dynamics exhibit 
universal scaling (§s 3,6,7)); and, relatedly, the behaviour of the renormalised level 
e* = €i + E^(0) (equation (2.14)). Unless stated to the contrary, all results given are 
representative of the full range < r < 1 considered in the paper. 

We begin by considering the GFL phase, U < Uc- The critical behaviour of 
the low-energy Kondo scale Wm = is illustrated in figure 3 for r = 0.3 and an 
asymmetry rj = 1 — 2\ei\/U = 0.2; the critical Uc{r, rj) 2± 12.2 being marked as a solid 
square. Specifically, {(jJm/Uy vs U (open circles) is shown, and seen to vanish linearly 
as U Uc—- The critical exponent a of §3 (w* ~ u" with u = \1 — U/Uc\) is thus 
a = 1/r, which agrees with the NRG calculations of [18] (the result is also derivable 
analytically for r <C 1, see §7.1). We find also that ui^ = Wm and the quasiparticle 
weight Z are indeed related by equation (3.8), namely that 



As argued in §3, the renormalised level ej* should likewise vanish as the transition 
is approached, from either phase. When considering the LM phase in particular it will 
in general prove helpful to define spin-dependent renormalised levels via 



by equation (2.11), from which the {e*^} and the renormalised level = Ci -I- E^(0) 
are related by 



5.1. Scales 





(5.21) 



Local quantum phase transition in the pseudogap Anderson model 16 



0.4 




Figure 3. Critical behaviour of the low-energy scale ui» in the vicinity of the 
QCP; illustrated for r = 0.3 and asymmetry = 0.2 (critical Uc — 12.2, solid 
square), ui^ vs U is shown: for U < Uc (GFL phase, uit, = ujni), [LOm/UY U is 
plotted (open circles); for U > Uc (LM phase, = uij^), |S^(0) - t,f(0)\ oc ojI 

is plotted (open squares). In either ease, ui* ~ |i7 — i7c|~ as the transition is 
approached. See text for further details. 



In the GFL phase (where symmetry is restored, S|'(0) = E|'(0)) the renormahsed 
levels all coincide, ej* = ef^. For the LM phase by contrast, e*,e*| and e^^ are in 
general distinct. 

Figure 4, again for r = 0.3 and rj = 0.2, shows the resultant [/-dependence of e* 
(open circles for U < Uc in the GFL phase, squares for U > Uc). The renormalised 
level is indeed seen to vanish as the transition is approached from cither phase, with 
an exponent ai = 1 (|e*| ^ agreement with the inferences drawn in §3. For 

the GFL phase in particular, and |e*| are thus related by equation (3.11), 

|e*| = kiul (and where the constant k = k(i]) is an odd function of the asymmetry ry, 
such that e* = by symmetry (§3) for the GFL phase of the particlc-holc symmetric 
PAIM, 7] = 0). For the LM phase U > Uc, figure 4 also shows the spin-dependent 
renormalised levels e*^ (triangles); each being seen to be sign definite and to vanish 
linearly as [/ Uc+- 

For the GFL phase, the low-energy (Kondo) scale cj* = arises as a direct 
consequence of symmetry restoration. The reader will however notice that the 
corresponding scale for the LM phase, in terms of which single-particle dynamics 
exhibit scaling and denoted by = wlj has not yet been specified. Here wc 
simply assert, and consider further in §s 6ff, that is related to the spin-dependent 
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Figure 4. Renormalised level e? = ej + S^(0) vs U in the vicinity of the QCP for 
r = 0.3 and r] = 0.2. GFL phase points {U < Uq) are shown as circles, LM phase 
as squares. In the LM phase the t / i-spin renormalised levels e*^ = e; + S^(0) 
are also shown, as up/down triangles. All vanish linearly as the transition is 
approached. 



renormalised levels by 

= WL oc lefi- - el^\-r = |Sf (0) - Sf (0)|^ (5.22) 

(or equivalcntly by cjl oc | e*^. | ; which equivalence follows since each e*^ vanishes 
linearly asU ^ Uc+, figure 4, and low-energy scales are equivalent up to an arbitrary 
constant. The specific proportionality employed in practice will also be specified in 
§6.2). Figure 3 also shows |S|^(0) — Sj^(0)| oc cj^ for U > Uc (open squares); the 
overall figure thus showing the i7-dependence of col in both phases. 

5.2. Phase diagrams 

We now consider representative phase diagrams arising from the LMA, obtained as 
described in §4. As usual the range < r < 1 is considered; and we show in particular 
that the essential characteristics of the phase diagrams divide into two classes, namely 
r < ^ and r > ^. 

Figure 5 shows behaviour typical of r < ^, the phase boundary in the {U, ei)-plane 
for a fixed r = 0.2. As ei is progressively increased, the critical Uc = Uc{r,ei) (solid 
line) first decreases, reaches a minimum at the particle- hole symmetric line Ci = —■^U 
(shown as a dashed line), before rising sharply at Cj ~ —5 and eventually turning 
back on itself. This has two interesting consequences. First that the LM phase is 
inaccessible by increasing U for any fixed ^ — 5 in the present example. Second, for 
fixed ei and due to the above 'turnback', there is at large U a re-entrant transition 
back to the GFL phase (itself seen most clearly in figure 18 below, where a low-r 
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Figure 5. LMA phase boundary in the {U, ei)-plane for a fixed r = 0.2 (solid 
line), representative of the general behaviour for r < i . The dashed line marks the 
paxticle-hole symmetric case: U = —26;. The corresponding MF phase boundary 
is shown as a dotted line. 



phase boundary determined analytically in §7.1 is shown). Both observations are in 
general accord with results of recent NRG calculations [36]; and indeed are seen even 
at the pure MF level of unrestricted Hartree-Fock, shown in figure 5 as a dotted line 
(although the LM phase is strongly exaggerated by pure MF, the many qualitative 
deficiencies of which are discussed e.g. in [15,25,27]). 

A rather more revealing portrayal of the phase boundary for fixed r is shown in 
figure 6 where, for r ~ 0.2, Uc = Uc{r,ri) is displayed as a function of the asymmetry 
T]. Specifically, 77 vs (8/7rr)J7^~^(r, 77) is shown; the constant S/nr simply serving as a 
convenient 'normalisation', it being known from previous LMA work [15] that as r ^ 0, 
U^~^{r,0) ~ 7rr/8 for the symmetric PAIM rj = (which result is exact [15,16]). In 
contrast to Uc{r, Cj) shown in figure 5, which is not single-valued in a (leading to 
re-entrance as above), figure 6 shows Uc{r,r]) to be a single- valued function of the 
asymmetry (which is found for all r); and that LM states arc always accessible upon 
increasing U for any fixed asymmetry |r;| < 1. In §7.1 the 77-dependence of Uc{r,rf} 
arising within the LMA will be discussed further and obtained analytically in the 
small-r limit (which result is also shown in figure 6). 

Figure 6 also shows clearly that particle-hole symmetry is optimal for LM 
formation — 77 = corresponding to the lowest critical Uc{r,r]). For the symmetric 
model it is however well known [10-20] that for r > 1/2 the GFL phase docs not 
survive for any non-zero U. It is for this reason that phase diagrams for r ^ 1/2 
are distinct, and figure 7 illustrates the point. For r = 0.7 the main figure (cf figure 
5) shows the critical Uc{r,ei) vs ei; from which in particular the re-entrant behaviour 
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Figure 6. Phase diagram in tlic (J7, r;)-planc for r = 0.2: t] vs {8/nr)U''~^ is 
sliown (solid line). Tlic corresponding r — + result, obtained analytically in §7.1, 
is also shown (dashed line); see §7.1 for full discussion. 
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Figure 7. Main figure: LMA phase boundary in the (C/, ei)-plaiie for r = 0.7, 
respresentative of the general behaviour for r > ^. The dashed line marks the 
symmetric limit, U = —26;; in this case (77 = 0) only LM states are accessible. 
Inset: phase boundary in the ({/, r))-plane, r] vs {S/'Kr)U^~^ being shown. 




Figure 8. r-dependence of the phase boundaries for fixed asymmetry rj: U~ 
vs r for (from top to bottom) rj = 0.1,0.2,0.3,0.4 and 0.5. See text for details. 
Since the wide-band Umit {D = oo) is employed here, the largest r considered in 
practice is r = 0.95. 



upon increasing U for fixed ej is again apparent. And the inset to figure 7 shows the 
?7-dependence of the phase boundary (in direct parallel to figure 6). The qualitative 
differences between the r < 1/2 and r > 1/2 phase boundaries are self-evident. 

The r-dcpendcncc of the phase boundaries is illustrated in figure 8, which shows 
U~^(r, rj) vs r for a range of asymmetries rj. For rj > the resultant phase boundaries 
have a maximum — in qualitative accord with the NRG results of [13] for [/ = oo 
— which moves to lower r with increasing r). For the symmetric case rj = the 
present approach is known [15] to be deficient as r ^ 1/2, in not capturing the slow 
(logarithmic) divergence of U~^{r — > |,0) seen in NRG calculations [12,13,16]; in turn 
reflecting the neglect of charge fluctuations in the present LMA (which are necessarily 
important in this limit, since Uc{r ^,0) 0). We do not however doubt at least 
the qualitative validity of the results shown for 77 > where Uc{r, rf) remains finite for 
all r considered. 

Finally, note from figure 8 that for small r the critical U~^{r,r]) oc r for fixed 
asymmetry, with a coefficient that is largest in the particle-hole symmetric case 
?7 = 0. We will revisit the matter analytically in §7.1, and show that this behaviour 
is intimately related to the exponential asymptotics of the ubiquitously non- vanishing 
Kondo scale characteristic of the r = metallic impurity model. 




Figure 9. GFL phase spectrum for r = 0.3, asymmetry r/ = 0.2 and interaction 
U = 11 (the critical Uc{r,ri) ~ 12.2); 7r£)(u)) vs w. Insets: spectral behaviour on 
the lowest energy scales. See text for discussion. 



6. Single- particle dynamics 

We turn now to LMA results arising for single-particle dynamics, on all relevant energy 
scales and including in particular the many different facets of their scaling behaviour 
as the underlying quantum phase transition is approached. The GFL and LM phases 
are considered separately in §s 6.1 and 6.2 respectively, and dynamics precisely at the 
QCP in §6.3. 

6.1. GFL phase 

A typical GFL phase spectrum is given in figure 9: TrD(oj) vs oj, shown for r = 0.3, 
an asymmetry r] = 1 + 2ei/U = 0.2 and an interaction strength U = 11 close to 
the critical Uc{r,ri) ~ 12.2. Note first that the change in number of electrons due to 
addition of the impurity is correctly riimp = 0, as required by the Freidel sum rule 
equation (4.8) (where e* = ei + S^(0) > for 77 > 0). By contrast, the local impurity 
charge n = 2 dLvD{uj) is n ~ 0.96, close to the Kondo limit whore n ^ 1; showing 
in turn (since rUmp = 0) that essentially one electron is transferred from the conduction 
band to the impurity upon its addition to the host. 

On high energy scales the expected Hubbard satellites arc seen in the spectrum, 
centred on a; ~ e; and Cj + C/; and correctly broadened two-fold over their counterparts 
at pure MF level for the physical reasons discussed in [15,25]. Spectral evolution with 
U is illustrated in figure 10 where, for the same r and 77, the spectra arc shown upon 
decreasing U further from the phase boundary into the GFL phase. With diminishing 
interaction the Hubbard satellites become progressively less pronounced as expected. 
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Figure 10. GFL phase spectral evolution with U, for r = 0.3 and rj = 0.2 
(the critical Uc ~ 12.2): -kD(lo) vs lo for U = 11,10,9,6 and 5 in obvious 
sequence. Inset: hD{lo) vs \lo\ on a logarithmic scale, for both a; > (solid 
line) and u) < (dashed), w = LOt (= LOm) is indicated by a vertical arrow. The 
function cos^(^r)|a;|~'' is also shown (dotted line); see text for full discussion. 



and ultimately absorbed into the central resonance structure. It is the latter that we 

now begin to consider in more detail, focusing first on the lowest energy behaviour. 

The left-hand inset of figure 9 shows the low-energy behaviour of the spectrum on 
a greatly expanded scale. It consists of two pronounced peaks of unequal intensity, one 
either side of the Fermi level lj = Q, with the more intense peak arising for > 0. This 
is quite general, applying for all r; > (where riimp = 0) and r, for U close to J7c(r, rf). 
For ?7 < this situation is reversed: under 77 —77 as discussed in §2, the spectrum is 
simply reflected, D(uj,ri) = D{—lo, —rj) (and rii^p = 0^2 — nimp = 2). The greater 
peak then lies below, and the lesser above, the Fermi level. For any asymmetry rj 0, 
where e* 7^ in the GFL phase, D{u)) has the ultimate low-w behaviour D{u)) oc |a;|''; 
given trivially (e.g. from equation (3.9)) by 

.Di.) ^ (6.1) 

and symmetric in u. This is illustrated in the right inset to figure 9 (and is of course 
in contrast to the symmetric PAIM [12,14-16,18,20] r? = (and ef = 0), for which 
D{w) oc \u)\~^ as |a;| — » 0, equation (3.5)). 

The qualitative features of the preceding paragraph are in agreement with results 
of recent NRG calculations [18]. We now show that they may be understood using 
the simple low-w quasiparticle form for the impurity spectrum, discussed in §3, 
that is characteristic of the GFL phase. Following that however, we point out the 
severe inadequacies of the quasiparticle form in understanding the GFL phase scaling 
spectrum on all relevant frequency scales, and in particular its essential irrelevance 
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in understanding the behaviour of the spectrum at the quantum critical point (which 
itself is considered in §'s 6.3,7.2). 

The quasiparticle spectrum is given in scaling form by equation (3.9), as a function 
of tD = w/w* with = Wm the low-energy Kondo spin-flip scale characteristic of 
the GFL phase (determined via symmetry restoration as described in §4); and with 
ki = u)l~'^/Z (equation (3.8)) a constant of order unity as the transition is approached, 
where and Z separately vanish (§5). The quasiparticle form is of course restricted 
by construction to <C 1 (which regime is now considered); and for present purposes 
equation (3.9) may be replaced with negligible error (as shown below) by 

nujlD{uj) — (6.2) 

[^-sgn(c.)/3(r)p|12 + p|2. 

where (3{r) = tan(^r) as usual (and with k = \el\/u!l likewise finite as C/ — > Uc{r, rj) — , 
as shown explicitly in §5). Equivalently and conveniently, equation (6.2) may be recast 
as 

'^l^'*!^^"^^ ^ [l-sgn(c.)/3(r)|.;'|1^ + kP'- ^^"^^ 

with (j' = w/lej"!^ (and where e* > has been taken, as appropriate to q > 
0). Equations (6.2) or (6.3) capture the low-energy 'two-peak' spectral structure 
illustrated in the insets to figure 9. Equation (6.1) is of course trivially recovered as 
\uj\ ^0; and from equation (6.3) the resultant spectrum has maxima at a; = ±T*, 
where 

T* = [|er|cos(fr)]^. (6.4) 

For any non-zero asymmetry 77, the low-energy scale T* is proportional to a;* and is 
thus equivalent to it (as follows from equation (3.11), noting that k (= k{r,r])) is odd 
in ?7 such that for the particle- hole symmetric model 77 = 0, = as in §3). We 
have here introduced T* simply to make connection to the NRG calculations of [18] 
for the asymmetric PAIM, where the low-energy scale was defined as the location of 
the spectral maxima. 

In figure 11, LMA results for Tr\e*\D{uj) as the transition is approached arc 
given, plotted vs oj/T* for a range of different r's (and noting that in all cases 
the range \u)\/T* < 10 shown corresponds to \w\/uJt, <C 1, since the constant 
k = \e*\/ujl ~ {T*/uj^Y is found in practice to be small compared to unity). We 
emphasise that the figure shows 'full' LMA results; that the spectra are of course 
scaling spectra, obtained as U ^ Uc{r,r])— and as such universally dependent on 
uj/T* (or equivalently uj) alone; and that the exponent b of §3 (equation (3.2)) is 
indeed b = r. As seen from figure 11 the asymmetrically disposed spectral maxima 
are indeed peaked at iv/T* = ±1, and for the particular case of r = 0.3 the figure 
compares explicitly the LMA results to the approximate quasiparticle form equation 
(6.3). This is not perceptibly distinguishable from the full LMA results, showing that 
equation (6.3) is indeed accurate (the same being true for the other r's shown). 

For sufficiently low oj/T* or ui, the GFL phase scaling spectrum is thus explicable 
in terms of the simple approximate quasiparticle form equations (6.2) or (6.3). An 
obvious question then arises: can this form be extrapolated beyond its strict confines, 
to infer the leading large uj/T* (or a)) form of the scaling spectrum which — as 
shown on general grounds in §3 (equations (3.2-4)) — determines the behaviour of the 
spectrum precisely at the quantum critical point where the low-energy scale oc T* = 




Figure 11. Low-frequency GFL scaling spectra: ■jT\e*\D(uj) vs lu/T* (with the 
low-energy scale T* defined in equation (6.4)), shown for r = 0.1 (dotted line), 
0.2 (short dash), 0.3 (solid), 0.4 (long dash) and 0.5 (point dash). For the case of 
r = 0.3, explicit comparison is also made to the approximate quasiparticle form 
equation (6.3); it is however indistinguishable from the full LMA results (which 
were obtained explicitly for r) = 0.2 but for »? ^ are in fact r^independent over 
the u;/T*-range shown, see equation (6.3)). 



0, and the generic form of which is given by equation (3.4)(7rZ?([/c; w) = C{r;r])\co\~'^)'! 
The reader will notice for example that the leading high frequency behaviour predicted 
from equations (6.2,3) is ttlvIDIu) ~ cos^(^r)|(2)|^'': from which the scale drops out 
of either side, resulting in an apparent QCP spectrum that is indeed of the required 
general form equation (3.4), but with a coefficient C{r,rj) =cos^(-|r). As seen from 
equation (3.5) the resultant QCP spectrum, if correct, would in fact coincide precisely 
with the non-interacting limit of the particle-hole symmetric (ry = 0) PAIM. 

The latter conclusions are not in fact correct, but the above discussion is 
nonetheless informative as now shown. Figure 12 gives the full LMA scaling spectrum 
for r = 0.3 and r] = 0.2: ttujID{lj) for w > (solid line) as a function of a) = w/w*, 
on a logarithmic scale encompassing the entire relevant Q range (and with a) = 1 seen 
as a small 'dip' in the spectrum, a well known and harmless artefact of the RPA- 
like form for IT"' employed here [15,25-27]). The figure also shows the approximate 
quasiparticle form equation (6.2) (dashed line), together with its asymptotic behaviour 
nu)lD{oj) ~cos^(^r)|ci;|~'" mentioned above (upper dotted line). The full LMA scaling 
spectrum is well captured by the quasiparticle form up to uj ~ 10~^, below which is 
indeed seen to exist an appreciable 'intermediate' a)-range over which the behaviour 
TrcjlD{u;) ~cos^(-|r)|£l;|~'' clearly arises. For larger frequencies however, the scaling 
spectrum clearly departs from '~cos^(^r)|d}|~'' behaviour; and for oj ^ 1 is seen to 




Figure 12. Full GFL phase scaling spectrum for r = 0.3 and rj = 0.2: ■kll)^D{ll)) 
vs u) = Lo/ujf {u)t = i^m) plotted for > on a logarithmic scale (solid 
line) and encompassing all relevant ui regimes. The approximate quasiparticle 
form equation (6.2) is shown as a dashed line; the upper dotted line shows the 
rj = non-interacting behaviour cos^(^r)|aj|~''. For |aj| 3> 1 the ultimate high- 
frequency behaviour of the scaling spectrum is also oc but with a different 
coefficient C{r,ri) that is r- and ?7-dependent; the lower dotted line shows the 
latter behaviour extrapolated back to lower lu. Full discussion in text. Inset: 
corresponding results for r = 0.7 and r] = 0.4, same notation as main figure. 



cross over into ultimate large-w behaviour that is again oc \uj\~^ — foUowing generally 
from the fact that the exponent 6 = r as explained in §3 — but with a prefactor 
C(r, ?/) that differs from cos^(^r) and is found in practice to be 77-dependent. 

The above behaviour is not particular to the case r = 0.3 illustrated, being found 
to arise throughout the range < r < 1 considered (the inset to figure 12 shows 
corresponding results for r — 0.7 and 77 = 0.4). The asymmetric GFL phase scaling 
spectrum thus contains two regimes of {Col^^ behaviour: that at lower uj (<C 1) 
corresponds to symmetric (rj = 0) non-interacting limit behaviour, while that arising 
for (D 3> 1 determines the QCP spectrum. This 'two \lu\^^ regime' behaviour is in fact 
known [20] to arise also in the pure symmetric PAIM, 77 = (where the cos^(^r)|a)|~'' 
regime persists all the way down to a; — 0, reflecting the exact low-w behaviour 
equation (3.5) arising in that case). Further consideration of the matter will be given 
in §6.3, and the (r, 77)-dependence of the QCP spectrum arising within the LMA will 
be determined analytically for small r in §7.2. 

Finally, we point out that figure 12 shows the scaling spectrum for w > 0, but 
that it is not of course symmetric in to (or uj /T*) as is obvious e.g. from figure 11. The 
sgn(a;)-dependence of the GFL phase spectrum is illustrated further in the inset to 
figure 10 above (for [r, 77] = [0.3,0.2] and U — 11, quite close to Uc — 12.2). Here the 
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Figure 13. LM phase spectral evolution with U, for r = 0.3 and r] = 0.2 
{Uc ~ 12.2); nD{uj) vs ui for U = 50, 20, 17, 14 (clockwise from top left). 



spectrum nD{u) vs \u)\ (i.e. unsealed) is again shown on a log-scale, for both w > 
(solid line) and u < (dashed). For each sgn(a;) the above 'intermediate' behaviour 
~cos^(^r)|cj|^'" arises (shown as a dotted line). Beyond that however, a crossover 
to ~ C{r,ri)\uj\^^ behaviour occurs; with the C{r,rj) dependent not only on r and rj 
but also for ?7 7^ on sgn(a)) (which behaviour we add is not precluded by the general 
considerations of §3, and which is likewise considered further in the following sections). 
Since the spectrum here is shown in unsealed form the non-universal Hubbard satellites 
are also ultimately visible on the figure, likewise asymmetrically disposed about w = 
because rj ^ 0. Wc also add here that while results from the present LMA are in 
general agreement with the NRG calculations of [18], there is one point on which we 
differ. That we find the coefficient C(r, 77) to depend on 77 (albeit fairly weakly), is 
consistent with a line of (ry-dependent) critical fixed points. Numerical RG calculations 
for the pseudogap Kondo model [13] appear by contrast to be consistent for any r] > 
with a single critical fixed point, and hence in effect no 77-dependence in C{r,ri). 



6.2. LM phase 

We turn now to the LM phase, considering first spectral evolution on 'all scales' as U 
is progressively decreased for fixed asymmetry 77, towards the LM/GFL transition at 
Uc{r,ri). This is illustrated in figure 13, again for r = 0.3 and r] = 0.2 (Uc — 12.2), 
although the behaviour shown is representative of the full r-range considered. TrD{uj) 
vs u! is shown, for U = 50, 20, 17 and 14. For the highest U shown, deep in the LM 
phase, the well developed Hubbard satellites dominate the spectrum entirely, although 
the ultimate low-w spectral behaviour vanishes as oc \w\'^ throughout the LM phase 
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Figure 14. LM phase spectrum nD{uj) vs w for r = 0.3, t] = 0.2 and interaetion 
U = 12.3 elose to the critical Uc — 12.2. Left inset: low-energy behaviour on an 
expanded scale. Right inset: nD(Lu) vs on a logarithmic scale, for both a; > 
(solid line) and u) < (dashed). Pull discussion in text. 



as noted in §3. With decreasing U however, a low-energy spectral structure is seen 
to develop in the vicinity of the Fermi level. It becomes increasingly pronounced as 
U Uc+ and as shown below is the precursor, in the LM phase, of the \uj\~'^ divergent 
behaviour characteristic of the QCP itself. 

For U = 12.3, close to Uc, ttD{uj) vs iv is shown in figure 14; the left inset showing 
the low-energy part of the spectrum. The latter is clearly redolent of its counterpart 
for the GFL phase (figure 9), the low-energy behaviour again consisting of two peaks 
of unequal intensity, the more intense peak arising for w > (and vice versa for 
Tj < 0). In contrast to the GFL phase however, the positions of the peak maxima 
are not symmetrically disposed about the Fermi level; reflecting the fact that in the 
LM phase, where symmetry is not restored, the spin-dependent renormalised levels 
e*^ = Ci + S^(cl; = 0) do not coincide (as shown e.g. in figure 4). In addition, and 
again in contrast to the GFL phase, we find within the LMA that while the lowcst-w 
spectral behaviour is oc as noted above, the coefficient thereof depends on sgn(a;). 
The origins of this are readily seen using equations (2.9,10), which give the leading 
low-w behaviour 

-kD{uj) ~ 2 r5n2 (6-5) 

with the LM phase S<t('^) given explicitly within the LMA considered here by equation 
(4.18). It is the first term of the latter that controls the low-w behaviour of ^^(a^) 
which, using equations (4.6,15), is found to be of form ^^(w) oc \ujy with a coefficient 
dependent on sgn(a;) (as shown explicitly for r? = in [15]); and which, via equation 
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Figure 15. Full LM phase scaling spectrum for r = 0.3 and rj = 0.2, obtained 
as (/ ^ Uc{r,T))+: itljID(lj) vs u) = lo/uj^ (a;* = wl) on a logajrithmic scale, for 
both a) > (solid line) and u) < (dashed). 



(6.5), generates corresponding behaviour in D{lo). 

The low-w peaks in figure 14 (left inset) mark the onset of a crossover from the 
asymptotic oc |a;|'' behaviour, to an \ui\~^ regime which persists essentiahy up to non- 
universal frequency scales on the order ~ 1 (= Aq). This is evident from the right 
inset to figure 14 where the spectrum is plotted on a logarithmic scale: the crossover 
from D{uj) oc \ujY to oc \u>\~'^ behaviour (likewise dependent on sgn(a;), albeit more 
weakly so) is clear, as too are the non-universal Hubbard satellites on the highest 
energy scales. 

The spectra shown in figures 13, 14 arc of course displayed on an 'absolute' scale, 
i.e. shown vs lo in units of the hybridization strength Aq = 1. But the LM phase is also 
characterised (§5) by a vanishing low-energy scale = Wl as the LM/GFL transition 
is approached, so the corresponding LM phase scaling spectrum is readily obtained 
in direct parallel to that for the GFL phase: on decreasing U progressively towards 
?7c-l-, steadily decreases, but for fixed r] the spectra again exhibit universal scaling 
in terms of a) = lo/loi, alone (and with non- universal features such as the Hubbard 
satellites thus 'projected out' as usual). As anticipated throughout we find that it 
is indeed 7r(jj££'(a>) vs u> which exhibits scaling (whence, as for the GFL phase, the 
exponent /; of §3 is again found to be 6 = r as required on general grounds); and with 
Wl given by the form equation (5.22). Specifically we take wl = to be given by 

= If [6*1 - 411" = If (0) - Sf (0)]|^ (6.6) 

where our (free) choice of prefactor, f , is chosen to make ready connection to the 
analytic results for small-r obtained in §7. The resulant LM phase scaling spectrum 
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for r = 0.3,7] = 0.2 is shown in figure 15, showing clearly both the ~ behaviour 
at low-(D and the crossover to the ultimate large-w behaviour oc 

6.3. Quantum critical point 

Which brings us to the quantum critical point itself. From the previous discussion we 
would naturally expect GFL and LM phase spectra arbitrarily close to the QCP to be 
indistinguishable for scales |a;| ^ w*. This is evident in figure 16 which shows spectra 
7rZ)(w) vs u) close to and on either side of the transition \Uc{'i' = 0.3, rj = 0.2) ~ 12.2]; 
for U = 12.0 in the GFL phase and U = 12.4 for the LM phase. On higher energy 
scales the spectra are indeed almost indistinguishable — limited only by the practical 
numerical inconvenience of dealing with the arbitrarily small but strictly non- vanishing 
scales in the immediate vicinity of the transition. 
The inset to figure 16 shows the corresponding a; > spectra on a logarithmic 
scale. The leading low-w behaviour ttD{uj) oc \uj\^ is clear, arising for each phase but 
with different coefficients for the; GFL and LM phases (and each such being r- and 
jy-dependent) . The figure also shows the 'kD{u) oc \co\~'^ behaviour arising for cj ^ u)^ 
in either phase; and, most importantly, that the coefficient of this |a;|~'" behaviour is 
identical for each phase — as argued on general gounds in §3. 

As /7 — » Uc{r, r])± and the low-energy scale — > 0, this power law persists 

to lower and lower frequencies, and alone survives at the QCP where = identically 
and the low-energy physics is scale-free. The LMA enables a direct determination 
of the spectrum precisely at the QCP U = Uc{r,r]) (numerically in general, and 
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analytically for small r as in §7). The resultant QCP spectrum for r = 0.3, = 0.2 is 

also shown in figure; from which the characteristic behaviour ttD{Uc', w) = C(r, r])\Lo\~^ 
as |a;| — > is clearly seen (and which in fact persists essentially up to non- universal 
scales u) ~ 0{1)). 

One obvious but important point should also be noted from the above discussion: 
close to the transition, and in either phase, acts as the natural crossover frequency 
scale to the emergence of common QCP behaviour for ::g> in single-particle 
dynamics, while for lower energies the spectral characteristics arc particular to the 
phase considered. This is of course the dynamical analogue, at T = 0, of the common 
schematic pertaining to static/thermodynamic properties at finite-T (see e.g. [21] or 
figure 1 of [19]), showing the crossover to QCP behaviour upon increasing T above 

6.4- Scaling: a closer look 

In the preceding sections we have focussed primarily on the asymmetric PAIM, rj / 0. 
The purpose of this section is three-fold. First, to compare results for ^ to those 
arising in the particle-hole symmetric limit r/ = [15,16,20]. There are of course 
important differences between the two — notably that for 77 = the LM phase alone 
arises for all r > ^ and any U ^ [10-20] — but, that aside, these are not as radical as 
might perhaps be inferred from the literature and the two may be readily understood 
on a common footing. 

In comparing dynamics for ry 7^ and 77 = 0, an obvious but important point 
should first be noted: that we have considered the scaling spectra in terms of cD = w/w* 
and that the low-energy scale is non-vanishing throughout both phases, including 
the particle-hole symmetric limit rj = [15,20]. In the GFL phase for example one 
might have chosen, as in §6.1, to express scaling in terms of w/T*; with T* cx |e*|~ 
given by equation (6.4) such that the low-energy maxima in the GFL phase spectra lie 
at ui/T* = ±1. This is the procedure adopted in [18] and — in general — is entirely 
equivalent to (w/w,) -scaling because T* cx w*. But the proportionality between the 
two, given (see equation (3.11)) from k{r,ri) = ^ l^)"" '^^ V: such that e* 

vanishes throughout the GFL phase in the symmetric limit, while Wh. = '-^m remains 
non-zero. For this reason, scaling in terms of is naturally required to consider the 
symmetric and asymmetric PAIMs on an equivalent footing. 

In comparing the GFL phase scaling spectra for different 77, it is also convenient to 
consider the so-called modified spectral function [14-16,20] J^{u)) =sec'^{^r)TT\ui\'^ D{uj), 
i.e. 

J^{uj) = sec^{^)r)\u}\^ XTTulD{co) (6.7) 

(likewise universally dependent on Q). For 77 = the resultant spectrum is 'pinned' 

at the Fermi level for all U in the GFL phase [14-16], J^{lu = 0) = 1 (as follows from 
equation (3.5)); while for any ^ 0, oc jw)^'' as \ui\ 0. The scale on which 

this distinction between the r? = and r? ^ scaling spectra arises is T* (oc w*) as 
discussed in §6.1 (figure 11). But in practice the ratio T*/uj^, is found within the 
LMA to be small (r*/w, ~ [k{r,r]]~ with k itself small compared to unity), so the 
qualitative differences between r] = and ^ are apparent only on the lowest-cD scales. 
This is seen clearly in the left inset to figure 17a which, for r = 0.3 and asymmetries 
7] = 0,0.2 and 0.5, shows J^{u)) vs Co for \uj\ <C 1. 
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Figure 17. Scaling spectra in both phases for fixed r (= 0.3) and varying 
asymmetry rj. A (upper): GFL phase scahng spectrum ^{uj) =scc'^ {^r)\uj\'^ X 
TrLj'^D(Lij) vs u) = u)/u)t (t<J* = LOn\) for asymmetries r} = (dotted line), 0.2 
(bold solid) and 0.5 (thin solid line). The main figure shows the central 'Kondo 
resonance'. The right inset shows the spectral 'tails' out to \Cj\ = 150, = 
(dotted line), i) = 0.2 (bold solid line for uj > bold dashed line for tj < 0) 
and ri = 0.5 (thin solid [dashed] line for w > [w < 0]). Left inset: scaling 
spectra for ;/ = 0,0.2 and 0.5 on the lowest scales (c/ figure 11). B (lower): LM 
phase scaling spectrum ixll>^D{ui) vs uj = uj/uit, (wt, = lui^) for rj = (dotted line), 
rj = 0.2 (bold solid line) and rj = 0.5 (thin solid line). Inset: spectra shown out 
to |a;| = 10^ for both tZ) > and < 0, same labelling as in figure 17A. 
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By contrast however, the strong similarities between particle- hole symmetric and 
asymmetric spectra in the GFL phase arc also evident in figure 17a. The main figure 
gives !F{uj) on the |a;| ~ 0{1) scale, showing clearly the low-energy 'Kondo resonance' 
in J^{uj) [14-16] (which becomes precisely the usual Kondo resonance of the metallic 
AIM in the r — limit [14-16]). The effects of non-zero rj, while present, do not lead 
to qualitatively different behaviour either on the ui ~ 0{1) scales or in the spectral 
'tails', shown out to tD ~ C(IO^) in the right inset. These comments hold also for the 
LM phase scaling spectra, t:luI^D{lu) vs Co = uj /uji^ being shown in figure 17b (again for 
r = 0.3 and r] = 0,0.2 and 0.5). In this case 7ra;££'(a;) oc as tD — > for both the 
symmetric and asymmetric cases. The spectra for different r} are thus qualitatively 
similar for all a), albeit that the extent of spectral asymmetry naturally increases with 

•n- 

In previous sections we have emphasised in particular the crossover from the 
oc \ujY behaviour of 'kuj^^D{u}) arising at low-a) to the ultimate largc-cj behaviour 
oc \(b\~^ that in turn controls the QCP spectrum. Our second point is simply to note 
that for the GFL phase, the asymptotic approach to the ultimate \Cb\~'^ form — where 
the corresponding J-^{u}) oc x loI^D{u}) thus plateaus to a constant — is quite subtle. 
In figure 17a for example, where the spectral tails of CF{uj) are shown out to \uj\ = 150, 
has not yet reached its plateau value but is still decreasing, albeit slowly; indeed 
from figure 12 it is clear the plateau docs not occur until cD ^ 0(10'^) for the case 
of r = 0.3 shown. The approach to the ultimate large-w behaviour of the scaling 
spectrum is thus slow, and in fact becomes progessively more so with diminishing 
r ^ 0. The nature of this behaviour (which is relevant to the GFL phase alone) 
is known from our recent work on the symmetric PAIM [20] and is modified only 
quantitatively for r] ^ 0: as r ^ 0, the ultimate large- |a)| behaviour arises only for 
\u}\ ^ e~ . For 1 <C <C e~ by contrast the decay of T{uj) towards its plateau value 
is logarithmically slow. This 'intermediate' large-tD regime progressively dominates 
the GFL scaling spectrum out to ever larger-tD scales as r — » 0, and for r = recovers 
precisely the known [26] logarithmic large-J; tails characteristic of the metallic AIM 
(for which -D(w) oc l/ln'^{\oj\) for \uj\ 1). While it is the small-r behaviour to which 
we turn in the final section, the reader is referred to [20] for further details on this 
point. 

As noted above, the ultimate large-w behaviour u)ID{lo) oc \oj\~'^ is not in fact 
reached in the a)-range shown in figure 17a for the GFL phase. In consequence none 
of the behaviour illustrated in the figure survives at the QCP itself. Our final remark 
here is simply to emphasise this point, evident also in figure 16 (inset): the Kondo 
resonance in !F{u)), illustrated clearly figure 17a and characteristic of the GFL phase, 
'collapses' at the QCP where its weight vanishes; as indeed remarked in [19], and 
which the above figures substantiate. 

7. Small r asymptotics 

Our aim now is to obtain analytical results for the critical behaviour of the low- 
energy scales u;* , the; phase; boundaries and single-particle scaling spectra, in the strong 
coupling regime where Uc{r, 77) » 1. This means in effect that we consider here small 
r <C 1, since oc r as r — > (see e.g. figure 8). 

In strong coupling, where the impurity charge n ^ 1, the low-energy behaviour 
of the pseudogap AIM maps under a Schrieffer- Wolff transformation onto the 
corresponding Kondo model [2], with host density of states /)host(t^) = Po\'^\^0{D—\u}\). 
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The exchange coupHng matrix element J = 2F^[l/|ei| + 1/{U — |ei|)], i.e. J = 
8V'^ /[U{1 — ry^)] expressed in terms of U and 77 = 1 — 2|ei|/C/; and the corresponding 
potential scattering strength is K = rjJ. Since the PAIM hybridization Ai(a;) = 
7rl^^Phost(i^) (with Ai(a;) given by equation (2.4)), it follows that 

1 

ttU 1 ~r;2 

from which strong coupling PAIM results can be transcribed directly to those 
applicable to the Kondo model. 

As a straightforward generalisation of [15,20] we first focus on {ECT(a;)} in the 
GFL phase (equation (4.14)), where in strong coupling Scr(w) takes the form 

^Au)^U^gZ{uj + <jLU^). (7.2) 

This follows from the fact that the spectral weight of Imll^ (ui) in strong coupling 
is confined asymptotically to w > with /Q°°(dw/7r)Imn+~(ci;) = 1 (which physically 
reflects saturation of the local moment, 1); and that its resonance is centred 

on the low-energy spin-flip scale on which scales the {G^{ijJ)} are slowly varying. 
From equation (4.14), equation (7.2) thus results. Equation (7.2) with uim = also 
holds for the LM phase, with uj^^ = reflecting as ever the local spin-degeneracy 
inherent to the LM state; as follows from equation (4.18) noting that [15,20] the 
LM poleweight Q — > 1 in strong coupling (and that the ^So-(w) contribution is both 
negligible in intensity and subdominant in frequency [15]). 

Reg^{u!) is given by the one-sided Hilbert transform equation (4.15), and its 
leading low-u behaviour in strong coupling (and hence for r — > 0) is readily obtained 
using the methods of [15,20,27]. The result is 



TTr 



1_ H 

A'' 



O{\oj\) (7.3) 



where A —imn{D,U) arises as a high-energy cutoff [15,20,27] and = 9a{e[) = 
[1 -I- cre|]~2, with e[ = ^e-jU and Ci determined as specified in §s 4.1.1 and 4.3. 
Equations (7.2) and (7.3) for T,^{uj = 0) wiU now be employed to determine the 
critical behaviour of = Wm, and consequent phase boundaries. We also note that 
the introduction of A =min(£', U) means the following analysis encompasses not only 
D ^ U {as exemplified by the wide band limit D = 00 employed explicitly in the full 
calculations of previous sections), but in addition the case U ^ D. As will be seen, the 
essential differences between the two naturally reside only in the precise dependence of 
the low-energy scales and phase boundaries [15,16] on bare model parameters; while for 
obvious physical reasons the scaling spectra for both phases are entirely independent 
of whether D ^ U ov vice versa. 



7.1. Statics: scales and phase boundaries 

We consider first the approach to the transition from the GFL phase, determining 
thereby the critical Uc{r,rf) for small-r and the critical behaviour of the Kondo spin- 
flip scale cjm = uj* as U ^ Uc — . A brief discussion is then given of resultant phase 
boundaries, and the corresponding low-energy scale = characteristic of the LM 
phase. 



Local quantum phase transition in the pseudogap Anderson model 



34 



The symmetry restoration condition equation (4.5) is a central element of our 

approach to the GFL phase. Using equation (4.9) it is given in strong coupling (where 
1) by S? (0) - Sf (0) = U; and hence from equations (7.2,3) by 



AX- 



Prom this the critical Uc{r,r]), where Wm — 0, follows directly from 

(X /Uc{r,r])Y^'^{'Kr/4:){J2a9'r)~^j whence vanishes linearly in r as r — > as 
noted in the discussion of figure 8 (§5.1). In the particle-hole symmetric (r/ = 0) limit, 
e( = (for all r and U, see §4) and hence g+ = 1 = g-; the result obtained previously 
in [15] is then recovered, U~^{r,r] = 0) = irr/S as r — > 0, which is exact [15,16]. More 
generally a knowledge of the as r ^ is required. But since U~^(r,r]) cx r as 
above, the e[ and hence go- may be replaced asymptotically by their corresponding 
r = values [27]; these are [/-independent and functions solely of the asymmetry r], 
given explicitly within the present LMA by [27] 

'^ = yB fiv') = h[^ + {l-v')-]. (7.5) 

(This result may also be obtained by detailed consideration of the limit r — > itself, 

analysis of the renormalised level e* = ei-l-S^(O) showing both that the ga given above 
hold as U ^ Uc{r,ri) and that the exponent for the vanishing of |e*| ~ u""^ (equation 
(3.10)) is indeed ai = 1 as in figure 4.) From equations (7.4,5) the critical Uc{r,r]) is 
thus given from 

A'' r^o 7rr (1 - ry^) 



(7.6) 



or equivalently 



mPoJ. ^ (7.7) 

for the critical exchange coupling of the corresponding Kondo model. 

We shall return to this below, but note first that the U Uc{r,r])— critical 
behaviour of the Kondo scale Wm follows in turn from equation (7.4) as 



The exponent for the vanishing of the low-energy GFL scale Wm = is thus 1/r; as 
found in figure 3 where, with A =min(Z), U) = U as appropriate to the wide band limit 
D = oo there employed, {u)*/Uy vs U was shown to vanish linearly in m = {1 — U /Uc) 
as u — > 0+. We add moreover that using equation (7.6) in equation (7.8) and taking 
the limit r — > recovers precisely the Kondo scale hitherto found in [27] for the r = 
metalUc AIM, 



A --"v 8 /(„») ; - "V p„jfm 

(with J = J{U, r]^) given by equation (7.1)). While recovery of an exponentially small 
Kondo scale from an approximate theory is non-trivial, we note that the exponent 
therein differs in general from the exact result [2] for the r = Kondo model by the 
factor of f{rf) G [1, ^]. It is thus as such exact only for the symmetric model where 
f{T] = 0) = 1; albeit that the correction is modest, / being slowly varying in ry and lying 
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Figure 18. Phase boundary in the (?7, ei)-planc for fixed r = 0.05 (solid line), 
obtained from equation (7.6) (with \ = U as appropriate to the wide band case). 
The corresponding result for r = 0.01 is shown for comparison (dashed line). 
Inset: r = 0.05 phase boundary for large negative e;, which is seen to approach 
asymptotically Uc = |ei|+- 



e.g. within 10% of unity for r] < 0.6. Since equation (7.6) (or equivalently equation 
(7.7)) likewise reflects the r — exponent, the 77-dependcncc of the corresponding 
exact result as r — » is obtained by replacing /(?7^) by 1 therein. 

The r ^ phase boundary arising from equation (7.6) is shown in figure 6 of 
§5.2 (with \ = U as appropriate to the wide-band limit), where it is compared to 
full numerical results for r = 0.2 and displayed as rj vs {8/nr)U^~^{r,r]) — which we 
note is equivalently limr^o[Uc{r,r] = 0)/Uc{r,'r])]. Here we add in passing that the 
resultant r ^ phase boundary between non-magnetic (GFL) and magnetic (LM) 
states, is redolent of the phase diagram obtained in Anderson's original paper on the 
r = metallic AIM (figure 4 of [1]). There is of course no transition for r = 0: the 
'phase diagram' of [1] was simply the boundary to local moment formation at MF 
level, the qualitative deficiencies of which are long since known; in particular that a 
spurious transition arises at a critical U^^{r = 0,r]) (with [7^^(0,0) = tt). But the 
similarities between the two arc nonetheless striking and the r ~ MF results for 
ri vs [U^^{0,r] = 0)/C/^^(0, r;)] ([1], figure 4), while wrong for r = itself, appear 
ironically to be qualitatively correct as r — > 0. 

Figure 18 by contrast shows the phase diagram in the (C/, ei)-planc arising from 
equation (7.6) (with A = U), for both r = 0.05 (solid line) and r = 0.01; the former 
showing again the re-entrant behaviour discussed in §5.2 in connection with figure 
5. The increasing stability of the GFL phase with diminishing r is clearly evident, 
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0.2 0.4 0.6 0.8 

r 



Figure 19. Full LMA phase boundaries, siiown as F(r,r)) = (8/7r)[/(»?^ )/(l — 
ri'^)]U^~^ vs r, for different asymmetries rj = 0.1, 0.2, 0.3, 0.4 and 0.5 (from top 
to bottom in sequence). As r —* 0, all collapse to the common form F{r,T]) ~ r 
(itself shown, dashed line). 



and expected physically since the Fermi liquid phase alone persists throtighout the 
entire {U, ei)-plane in the r = limit of the metallic AIM. The inset to the fignre also 
shows the phase boundary (for r = 0.05) out to large negative ei, from which it is seen 
that the critical Uc asymptotically approaches |ei|+ (shown as a dotted line) — again 
physically natural, this being the condition for local moment formation in the atomic 
limit. 

In figure 19 the r-dependence of full LMA phase boundaries is shown for a range 
of asymmetries ri, and the wide-band limit employed in previous sections. Specifically, 
F{r,r]) = (8/7r)[/(?7^)/(l — r]'^)]UJ^~^{r,T]) vs r is shown, which from equation (7.6) 
should behave as F{r,r]) ~ r as r ^ such that phase boundaries for difi'erent r 
then collapse to a common curve. This behaviour is indeed seen in figure 19 and 
we note that, in line with results from a previous NRG study [13], the resultant 
phase boundaries have little asymmetry dependence up to r ~ 0.25. The asymptotic 
behaviour F{r,r]) ~ r is also shown on the figure and seen to be reached in practice 
below r ^ 0.1 or so (which range we add is known from a previous LMA and NRG 
study [16] of the symmetric limit to be enhanced further if one works instead with a 
finite bandwith D, with U/D > 1). 

We now consider briefiy the corresponding U Uc+ critical behaviour of the 
low-energy scale = characteristic of the LM phase, given by equation (6.6) 
(the resultant phase boundary Uc (r, ij) is of course correctly independent of the phase 
from which it is approached) . This is straightforward, since in strong coupling (where 
ImI ^ 1), tf{0) - tf{0) = (0) - T,f{0) - U from equation (4.9); and S^(0) 
is given from equations (7.2,3) with u>m = throughout the LM phase, whence 
Sf (0) - Ef (0) = {Ay/nr) J^a 9<t - U. But from equation (7.4) {AY /-Kr) 9<t = Uc, 

so (0) - Sf (0) Uc(U/Uc - 1) as t/ ^ Uc+, with C/c oc l/r as r ^ from 

equation (7.6). Hence wl oc (r|S?'(0) — E|'(0)|)^/'' given by equation (6.6) indeed has 
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the critical behaviour 

with exponent 1/r, as shown numerically in figure 3 (§5.1). 
7.2. Scaling spectra and the QCP 

We turn finally to the scaling spectra characteristic of the GFL and LM phases (figures 

12, 15), in particular the leading large uj = uj/uj^, behaviour thereof that determines 
the spectrum at the quantum critical point itself (§6.3). 

The impurity Green function is given by equations (2.9,10), which may be 
expressed quite generally as G{uj) = IX^ctI^^ ^ ^{^) ^ ~ (^^(w) — T,c{0))]~^, 
where we use Scr(u;) — S(j(0) = So-(w) — Scr(O) (trivially, from equation (4.9)). Here 
as usual, e*^ = ei + 2^(0) is the renormalised level (equation (5.20)); with the e*^ = ej* 
independent of spin cr for the GFL phase where symmetry is restored (equation (4.5)), 
but ^ e* for the LM phase, as shown in §5.1. As discussed in the preceding sections, 
(i) the scaling spectra for the GFL or LM phases are obtained by considering finite 
w = w/w* (with (J* = Wm or wl as appropriate), in the formal limit 0; and (ii) it is 

u)lG{ijj) (as opposed to G{u!) itself) which exhibits universal scaling in terms of uj. This 
follows directly tolGico) ^ i EM""^^ - <^7''M^) - ^^{^1 + (S.(co) - S,(0))]-i, 
in which a'^~''|ci)| may be neglected entirely for < r < 1 as considered. By contrast, 
from equations (2.4,5) for the hybridization, uj~'^A{uj) is non-vanishing in the scaling 
regime and given by w^^A(w) = — sgn(a;)[/?(r)-|-i]|tD|'" (independently of the bandwith 
D). The scaling spectrum 'jtu)ID{u) = —sgn{u)wllmG{u)) is thus obtained from 

-1 



sgn(o.)[/3(r) +i]|c:.r " ^ " <'{^Au^) " S.(0)) 



(7.11) 



which we add holds generally for all < r < 1 considered; and where e^^/ujl as 
^ has been shown in §5.1 to tend to a non-zero constant for both the GFL and 
LM phases (denoted therein as k{r,ri) for the GFL phase). 

Since ulGioj) scales, it follows from equation (7.11) that uj^^'lYiaiuj) — Scr(O)] 
is likewise universally dependent solely on uj — lu/lo,:. In the strong coupling, r — > 
regime of interest here, this may be shown directly using equations (7.2,3) (together for 
Ej,(w) with Img^iuj) = -sgn{uj)TrD"jLj)e{±uj) with D°{uj) the MF spectral density). 
For the GFL phase one thereby finds 

uri^^iu) - E^(0)] = -a^g, [\co + - 1] (7.12a) 

w-'-S^(w) = 4g, e{-[l + aoj])\u> + a]'' (7.12b) 

(noting that I]J^(0) = 0); where the = gcr{r]) are given explicitly by equation (7.5) 
and = Wni- For the LM phase by contrast, 

coZ^i^^iu;) - S^(0)] = -a—g^ (7.13a) 

LJ-^Ki^) = 45a e(-aw)|a;r (7.13b) 

with Lu^ = t^L. Equations (7.12,13), which we emphasise hold for all uj (and r ^ 0), are 
indeed of the required universal form; and for 77 = 0, where the ga = I, reduce to the 
results of [20] for the symmetric PAIM. Together with equation (7.11) they determine 
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the full (D-dependence of the scaling spectra for the respective phases. The resultant 

range of behaviour is rich, encompassing for example the intermediate logarithmic 
regime in the GFL phase mentioned in §6.4, and may again be obtained following the 
arguments of [20]. 

Here we consider solely the ultimate large-(D behaviour of the scaling spectra for 
either phase, relevant to the QCP itself. This arises for \u!\^ 1, and for r ^ (where 
/3(r) =tan(|r) ^ 0) equations (7.11-13) yield the requisite behaviour 



given explicitly for positive uj (for < 0, g+ and g_ arc simply interchanged). The 
result equation (7.14) applies to both the GFL and LM phases; and cancels out 
of each side of the equation to give the behaviour 7tD{u>) = C{r,r])\Lj\~'^ precisely at 
the QCP (which in that case extends down to uj = (§6.3)). As foimd numerically 
the QCP coefficient C{r,r]) is indeed both r- and Ty-dependent; with C(r,r]) ~ as 
r — > and as such from equation (7.6) proportional to l/[J7c(r, rj)]^ (or equivalently 
[PqJc^)i in which the interacting nature of the QCP is directly manifest. 

8. Summary 

We have developed in this paper a non-perturbative local moment approach to the 
T = pseudogap Anderson impurity model, with a conduction electron density of 
states proportional to \lo\^ and itself a longstanding paradigm [6-20] for understanding 
the topical and important issue of local quantum phase transitions. The transition 
here is between a generalised fermi liquid phase, perturbatively connected to the non- 
interacting limit, in which the impurity spectrum exhibits a Kondo-like resonance 
indicative of underlying local singlet behaviour; and a doubly degenerate local 
moment phase in which the impurity spin remains unquenched. In addition to phase 
boundaries, and associated critical behaviour of relevant low-energy scales, local single- 
particle dynamics in the immediate vicinity of the transition have been considered: in 
the GFL and LM phases separately, and at the quantum critical point itself where the 
Kondo resonance has just collapsed and the problem is scale-free. Notwithstanding 
the range of vanishing low-energy scales as the transition is approached, and their 
importance in understanding local spectra, dynamics close to the QCP in either phase 
arc nonetheless controlled by a single low-energy scale cj, and as such scale universally 
in terms of w/w* . A determination and understanding of these scaling dynamics on 
all w/cj* scales, including their continuous evolution 'from the tails' to pure QCP 
behaviour, has been a primary focus of the work. Approximate though it is the LMA 
is to our knowledge one of the few theoretical approaches currently capable of handling 
these issues; and the description arising from it is both rich and arguably a good deal 
more subtle than appears to have been appreciated hitherto. Much of course remains 
to be understood, not least the question of dynamics at finite-temperature, aspects of 
which will be considered in a forthcoming paper. 
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Local quantum phase transition in the pseudogap 
Anderson model: scales, scaling and quantum 
critical dynamics 

Matthew T Glossop and David E Logan 

Oxford University, Physical and Theoretical Chemistry Laboratory, South Parks 
Road, Oxford OXl 3QZ, UK 

Abstract. The pseudogap Anderson impurity model provides a paradigm 
for understanding local quantum phase transitions, in this case between 
generalised fermi liquid and degenerate local moment phases. Here we develop 
a non-perturbative local moment approach to the generic asymmetric model, 
encompassing all energy scales and interaction strengths and leading thereby to 
a rich description of the problem. We investigate in particular underlying phase 
boundaries, the critical behaviour of relevant low-energy scales, and single-particle 
dynamics embodied in the local spectrum D{u>). Particular attention is given to 
the resultant universal scaling behaviour of dynamics close to the transition in 
both the GFL and LM phases, the scale-free physics characteristic of the quantum 
critical point itself, and the relation between the two. 



Local quantum phase transition in the pseudogap Anderson model 



2 



1. Introduction 

The essential physics of a single magnetic impurity coupled to a fermionic host is 
embodied at its simplest in the Anderson impurity model (AIM) [1]: a correlated, 
non-degenerate impurity with local interaction U, hybridizing to a non-interacting 
host band with density of states /9host('^) (for a comprehensive review see [2]). Yet 
its simplicity is nominal — even for the conventional case of a metallic host, the 
basic model for imdcrstanding magnetic impurities in metals [1,2] and highly topical 
again in the context of quantum dots for example [3,4], or surface atoms probed by 
scanning tunneling microscopy (STM) [5]. Here the essential strong coupling (large- 
U) behaviour is that of the Kondo effect, characterised by a single low-energy Kondo 
scale ujk, manifest vividly in the Kondo resonance in the local single-particle spectrum 
D{w), and rendered no less subtle for being by now quite well understood. 

There is however one sense in which the metallic AIM is 'simple': it lacks a 
non-trivial quantum phase transition. The impurity spin is quenched ubiquitously by 
coupling to the low-energy degrees of freedom of the metallic host, the system is a 
Fermi liquid for all interaction strengths U and the Kondo scale, while exponentially 
small in strong coupling, never quite vanishes. 

That situation changes markedly if the host contains instead a pseudogap in 
the vicinity of the Fermi level [6] (cu = 0), namely /Ohost(^) oc \lu\^ with r > 
(and r = corresponding to the normal metallic AIM). This pseudogap Anderson 
impurity model (PAIM) has been studied extensively in recent years [6-20], particularly 
for large-C/ where the low-energy subspace maps onto that of the associated Kondo 
model, and from which work much progress and insight has accrued. The resultant 
depletion of host states around the Fermi level leads in general to the occurrence [6- 
20] of a quantum phase transition at a finite, r-dcpendent Uc- Above the quantum 
critical point (QCP) a degenerate local moment (LM) phase arises, with a residual, 
unquenched impurity spin. For U < C/c by contrast, and perturbatively connected to 
the non-interacting limit (U = 0), resides a 'strong coupling' or generalised Fermi 
liquid (GFL) phase in which the impurity spin is locally quenched and a Kondo 
effect manifest, reflected in a low-energy Kondo scale ojk- In the vicinity of the 
transition each phase is in fact characterised by a low-energy scale [17-20], here denoted 
generically by lu^; corresponding to the Kondo scale lu^ = luk for the GFL phase and 
its counterpart = wl [20] for the LM phase. Most importantly, and symptomatic 
of the underlying transition, the a;*'s vanish as the QCP is approached. 

Although the PAIM serves in general as a paradigm for understanding local 
(boundary) quantum phase transitions [21], many examples of specific systems with 
a power-law density of states have also been discussed. These include various zero- 
gap semiconductors [22], one-dimensional interacting systems [23], and the important 
problem of impurity moments in d-wave superconductors [17,18], notably Zn- and Li- 
doped cuprate materials. Likewise it has recently been shown [19] that the pseudogap 
Kondo model exhibits critical local moment fluctuations, and an associated destruction 
of the Kondo effect, very similar to those present at the local QCP of the Kondo lattice 
[24] that has been argued to explain the behaviour of heavy fermion metals such as 

CeCu6-a;AUa;. 

While great progress has been made in understanding the PAIM [6-20], much 
nevertheless remains to be understood; particularly regarding the quantum critical 
behaviour of dynamical properties close to and on either side of the transition, and 
precisely at the QCP itself. Since the low-energy scale is arbitrarily small (but non- 
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zero) sufficiently close to the transition, the frequency {u>) dependence of dynamics 
should scale in terms of it — i.e. exhibit miivcrsal scaling as a function of uj/lu^ — 
with resultant scaling behaviour that will depend on the particular phase considered, 
GFL or LM. What precisely is that behaviour for the selected dynamical property? 
And how from it docs (or can) one imderstand the scale-free physics characteristic of 
the QCP itself, where = identically? These are tough issues, especially from an 
analytical perspective, and some key aspects of which have recently been uncovered 
via numerical rcnormalization group (NRG) calculations [16,18,19], the local moment 
approach [15,16,20] and, for the multichannel pseudogap Kondo model, a dynamic 
large- A'' approach [17]. Single-particle spectra embodied in the impurity D{u}) provide 
an obvious example of dynamics [12, 14-20]; and an important one given e.g. the direct 
manifestation in the GFL phase of low-energy Kondo physics embodied in the Kondo 
resonance, as well as their more general relevance to STM and related experimental 
techniques. 

In this paper we develop a local moment approach (LMA) [25-30,15,20] to the 
PAIM at T = 0, for the general case of the particle- hole asymmetric model (noting that 
the corresponding Kondo model thus contains both exchange and potential scattering) . 
Single-particle dynamics form a central focus of the theory albeit, importantly, that 
its intrinsic formulation also enables phase boundaries between the GFL and LM 
phases to be determined, as well as the critical behaviour of associated low-energy 
scales [15,20]. Although developed originally to describe metallic Anderson impurity 
models, at T = [25-27], finite temperatures [28] and magnetic fields [29], the approach 
is quite general; readily extended for example to encompass lattice-based models 
within the framework of dynamical mean-field theory, such as the periodic Anderson 
lattice model [30]. It has also been applied previously to and has led to a number of 
new predictions for, the symmetric limit of the PAIM [15,20]; rather successfully, as 
judged by subsequent comparison to NRG calculations [16]. The asymmetric PAIM 
considered here is however another matter: the symmetric model, recovered though 
it is as a particular limit of the present work, has long been argued [13,17,18] to be 
a somewhat special case, in some sense divorced from generic asymmetric behaviour 
(although our own caveats on the matter will become clear in due course). 

The LMA, while certainly approximate, is also well suited to the present 
problem. Its intrinsically non-perturbative nature enables it to handle strong electron 
correlations and the attendant, central low-energy scales; while at the same time 
permitting access to dynamics on all energy scales, whether it be on the lu^ scales 
relevant to universal scaling behaviour, or scales on the order of a; ~ i7 relevant to 
the high-energy Hubbard satellites. Importantly too, and in contrast to difficulties 
inherent to a number of other potential theoretical approaches, the method handles 
both the GFL and LM phases on an essentially equal footing (the LM phase, for 
example, is certainly not a trivially decoupled/isolated spin). 

The paper is organised as follows. After the requiste background, we consider in 
§2.1 an obvious question: since the GFL phase is pertubatively connected to the non- 
interacting limit, satisfying as such the Luttinger theorem [31], what is the analogue of 
the resultant Friedel sum rule? That is entirely familiar for the conventional metallic 
AIM (r = 0) [2,32], and its counterpart for the PAIM plays a central role in the 
present work. In §3 another important issue is considered: what can be deduced on 
general, minimalist grounds about the behaviour of the scaling spectra close to the 
quantum phase transition; in both the GFL and LM phases and at the QCP itself? 
Such considerations are revealing, and dictate constraints to be met by any credible 
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approximate theory. They show moreoever that a variety of low-energy scales in 
addition to a;* enter the scaling behaviour of dynamics the quasiparticlc weight Z 
and renormalised level ej* for the GFL phase, together with the latter's spin-dependent 
counterparts for the LM phase. Each vanishes as the transition is approached, with 
distinct but related critical exponents; but nonetheless enter dynamics in such a way 
that universal scaling remains determined solely by the key low-energy scale w*. 

The LMA itself is developed in §4, focusing on the notion of symmetry restoration 
that is central to the approach [15,20,25-30], together with what likewise proves 
important in dealing with the asymmetric PAIM, namely self-consistent satisfaction 
of the Priedel sum rule/Luttinger theorem. Results arising from the approach are 
presented in §s 5-7, beginning with the critical behaviour of relevant low-energy scales 
(§5.1) and a determination of representative phase diagrams for the problem (§5.2). 
Single-particle dynamics are considered in §6, on all energy scales whether universal or 
non-universal; but with a dominant focus on scaling spectra in the GFL and LM phases 
— on all lv/lo^ scales — as well as the scale- free dynamics arising at the QCP. Finally, 
we consider in §7 the regime of small r <C 1, obtaining explicit analytic results for the 
critical behaviour of low-energy scales, resultant phase boundaries and single-particle 
scaling spectra. The paper concludes with a brief summary /outlook. 

2. Background 

With the Fermi level taken as the energy origin, the Hamiltonian for the spin-| AIM 
[1,2] is given by 



The first term describes the non-interacting host with density of states Phost('^) = 
J2k ~ ^k), while the second refers to the correlated impurity with bare site-energy 
Ci and on-site interaction U. The third term is the one-electron host-impurity coupling 
embodied in Vifc, taken for convenience as Vn, = V (although that is not required). 
The particle-hole asymmetry of the model is specified by the asymmetry parameter 



the impurity charge n = '^^{fiia) = 1 for all U. In strong coupling, where the relevant 
low-energy sector of the AIM maps under a Schrieffer- Wolff transformation onto a 
Kondo model containing both exchange (J) and potential (K) scattering contributions 
[2], the asymmetry rj is equivalently the ratio of the potential/exchange scattering 
strengths: r] = 1 - 2\ei\/U = K/J (see e.g. [27]). 

We focus here on the impurity Green function G(a;)(<-> G{t) = — i(T (cio-(i)cj^))), 
and hence the local spectrum D{lj) — -'^sgn{uj)lm.G{(jj). In addition to intrinsic 
interest in single-particle dynamics per se, this also provides a direct route to both 
the location of phase boundaries between the GFL and LM states [15] and the critical 
behaviour of associated low-energy scales, as will be shown in §s 4,5. 

G{lo) is conventionally expressed as 




(2.1) 




■i = --f and 



(2.2) 



G{oj) = [lu+ -€i- A{uj) - E(w)] 



-1 



(2.3) 
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with u}~^ = u) + iO+sgn(a)) and S(a') = — isgn{u)T} (uj) the usual single 

sclf-cncrgy; the corresponding non-interacting Green function is denoted by gi^uj). 
All effects of host-impurity coupling at one-electron level are embodied in the 
hybridization function A{uj) = Y^^^iki^'^ ~ ^k]~^ (= Ar — isgn(a;)Ai(w)). For 
the PAIM considered, Ai{uj) = 7rF^/9host('^) is given explicitly by Ai{uj) — 
Ao{\uj\/Aoy0{D — with Aq the hybridization strength and D the bandwidth 
{0{x) is the unit step function). That this form is naturally simplified is of course 
immaterial to the relevant low-energy physics of the problem in the scaling regime of 
primary interest. We take Aq = 1 as the basic energy unit in this paper, i.e. 

Ai{lu) :^\uj\''e{D -\uj\) (2.4) 

with r > for the pseudogap AIM (and the limit r = corresponding to the normal 
metallic model). The real part Aji{u) follows from Hilbert transformation and, for 
later reference, has the low-o; behaviour (for all r > 0) 



Anioj) = -sgn(a;) J f3ir)\u;\^ + + O 

I 7r(r ~ 1) D 



D 



(2.5) 



with (3{r) = tan(^r) used throughout. Notice from this that the wide-band limit 
Z) ^ 00, so frequently employed for the metallic AIM (see eg [2]), can be taken for 
r < 1, whereupon AR(a;) = —sgn{uj)(3{r)\uj\'' for all lo. We consider < r < 1 in this 
paper, which range is assumed unless stated to the contrary; and for that reason will 
employ the wide-band limit. This is simply a matter of convenience, is again irrelevant 
to the key low-energy behaviour of the system and can be relaxed trivially. 

The impurity charge n is of course related to the local spectrum by n = 
2 J^^dw D{u)). The 'excess' charge nimp will also prove important to the present 
work, for reasons given in §2.1. This is defined as the change in the number of electrons 
of the entire system due to addition of the impurity to the host, and is given generally 

by 



J — <x 



du Ap{uj). (2.6) 



Here Ap{uj) (= p(a;) — phost{<^)) is the excess density of states, itself given by (see e.g. 
[2]) 



Ap{lj) = — — sgn(a;)Im \ G{lo) 
n 



dA{w) 



dui 



(2.7) 



such that riimp is related to n by 



ni^p = n- h^n J du G{u)^^^. (2.8) 

From the above, for specified gap index r, the PAIM is characterised by any two 
of the 'bare' parameters U, Cj and r] = 1 + 2ei/U (together if appropriate with the 
bandwidth D). In practice, as now explained, we choose to consider ry and one or other 
of U and ei as the basic parameters. For the r = metallic AIM we have recently 
considered [26,27] the scaling behaviour of -D(u;) arising in the strong coupling Kondo 
regime, where Ci < such that |ei| ^ 1 and U — \e\\ ^ 1. Here the Kondo scale 
is exponentially small. For fixed asymmetry 77 [27], D{u}) then exhibits scaling in 
terms of a) = w/wk- That is, on fixing rj = 1 — 2\ti\/U and progressively increasing 
the interaction U, D{u) is universally dependent on a) = w/wk alone, and otherwise 
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independent of the bare model parameters (which themselves determine the Kondo 
scale). The important point here is that this arises for fixed asymmetry t/ = K/J 
(and in this case only for fixed 77, i.e. universal scaling does not arise upon increasing 
U for fixed |ei| » 1 [27]). In consequence, a family of scaling spectra arises, one for 
each asymmetry 77. This is not widely appreciated (although in the related context of 
d.c. transport it was understood early by Wilson [33]), but that it must arise is nearly 
obvious. For if the scaling spectra for the generic asymmetric AIM were 7?-independent, 
they would coincide with that for the p-h symmetric model 77 = (= K/J) and as 
such be strictly symmetric in u> about the Fermi level w = 0; the implausibility of 
which is at least intuitive. 

For the pseudogap model considered here, the quantum phase transition between 
GFL and LM phases is naturally characterised by a low-energy scale — the Kondo 
scale wk for the GFL phase and its counterpart wl [20] (discussed in §6ff) for the 
LM phase that vanishes as U Uc{r)± and the transition is approached. In 
consequence the spectra likewise exhibit scaling in terms of d; = w/w* in the vicinity 
of the transition. As discussed in §6ff, we again find this scaling to be ry-dependent. 
For this reason, and to enable continuity as r ^ to the above scaling behaviour 
for the metallic AIM, we choose to regard ry and one or other of U and Ci as the 
basic model parameters. The minimal range of 77 that need be considered is also 
readily established, for under a general p-h transformation of the Hamiltonian it is 
straightforward to show that ei — [ei + U] and -D(a;; Ci) = D{—lo; — [ei -|- U]), hence 
n(ei) = 2 — n(— [ci + [/]) for the impurity charge (and likewise for the excess charge 
n\aii>)- From this it follows that only ei > —\U need be considered, corresponding to 
77 > and which range is assumed from now on. 

To capture the underlying phase transition intrinsic to the PAIM, an inherently 
non-perturbativc approach is obviously required. That poses severe and, to 
our knowledge, insurmountable problems for any approach based directly on the 
conventional single self-energy Y,{uj). But such an approach is not obligatory, S(w) 
merely being defined by the Dyson equation implicit in equation (2.3) for G{uj). For 
the doubly degenerate LM phase in fact, a two-self-energy (TSE) description would 
be physically far more natural (as is obvious by considering e.g. the atomic limit, an 
extreme case of the LM phase). Here the rotationally invariant G{w) is expressed as 



in terms of self-energies Scr(w) (= S^(a;) — isgn(a;)Sj,(a;)) that are in general distinct 
and spin-dependent. The TSE description is in general a necessity and not a luxury 
for the LM phase. It is precisely this framework that underlies the LMA [25-27] where, 
as explained further in §4, it is employed for both the LM and GFL phases and in 
consequence is able to capture simultaneously both phases and hence the transition 
between them. The conventional single self-energy can of course be obtained as a 
byproduct of the TSE description. It follows from direct comparison of equations 
(2.3,2.9,2.10) and is given for later use by 



G{u;) = \[G^{uj) + G^{uj)] 



(2.9) 



with the G„{uj) given by 



G,(w) = [u;+ - e; - Aiuj) - E,(a;)] 



-1 



(2.10) 



S(a;) 




i(|]T(^)-Ex(c.)) 



(2.11) 



with g{uj) = [uj 
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Before turning to the LMA itself we consider two important questions. What 

can be deduced on general grounds about the behaviour of the GFL phase? That is 
considered in §2.1. Second, as considered in §3, what may be inferred via general 
scaling arguments about the behaviour of the scaling spectra close to the phase 
transition; in both the GFL and LM phases and at the quantum critical point (QCP) 
itself, where the low-energy scale w* = identically? 



2.1. GFL phase: Friedel sum rule 

The GFL phase is perturbatively connected to the non-interacting limit, which in 
essence defines that phase. This adiabatic continuity requires that the Luttinger 
integral vanish [2,31], i.e. that 

doj G{u) = 0. (2.12) 



= lm f 

J —cx 



From this, using only S'(0) = 0, the Friedel sum rule follows directly. Given explicitly 

by 

nimp = 1 - -tan-i ( -r^ ) (2.13) 



it relates the excess charge riimp to the 'renormalised level' e*, itself defined by 

er = ei + SR(0). (2.14) 

For the r = metallic AIM, where Ai(0) = Aq (=1), this gives the Friedel sum rule 
in the famihar form e* = tan[|(l - nj^p)] [2]. For the PAIM by contrast, Ai(0) = 0; 
and in consequence nimp is restricted to the integer values 

'o - t^n (2-15) 
2 : ef < ^ ' 

according to whether the renormalised level lies above or below the Fermi level. 
This behaviour, which we reiterate is a consequence of adiabatic continuity and the 
Luttinger theorem, is physically natural: one expects intuitively that the change in 
number of electrons due to addition of the impurity should be integral. It is in fact 
the more familiar metallic AIM, where from equation (2.13) above riimp is in general 
non-integral, that is the exception; reflecting in that case the existence of a finite 
density of conduction band states arbitrarily close to the Fermi level. 

The above conclusions may be checked order-by-order in perturbation theory (FT) 
about the non-interacting limit (as we have confirmed explicitly at simple first and 
second order levels) : the resultant ni„ip calculated from spectral integration (equation 
(2.6)) indeed satisfies equation (2.15), and Jl = 0. However the Friedel sum rule and 
hence Luttinger theorem is not satsified by approximate theories in general, save for 
the p-h symmetric case where it is guaranteed by symmetry (and for which Ui^p = 1 
regardless of whether the GFL or LM phase is considered). For the metallic AIM 
that is well known to be the case for the non-crossing approximation (NCA) [34] , and 
is true even for second order PT employing Hartree renormalised propagators [35], 
usually but loosely referred to as second order PT: the resultant nimp determined by 
spectral integration does not in general concur with that implied by the Friedel sum 
rule (whether for the normal metallic model or the PAIM). It is however obviously 
desirable that the Friedel sum rule be satisfied if possible, and in §4.1 we incorporate 
it straightforwardly within the LMA. 
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3. Scales and scaling: general considerations 

Close to and on either side of the quantum phase transition intrinsic to the PAIM, the 
problem is characterised by a single low-energy scale denoted generically by u* that 
vanishes as the transition is approached, U — > Uc{r)± [17-20]. Given this we now use 
general scaling arguments to make a number of deductions about the scaling spectra, 
in both the GFL and LM phases and at the quantum critical point (QCP) itself 
(where oj^ = 0). We take it for granted in the following that the QCP is continuously 
accessible from either the GFL or LM phases as ^0. Although natural, implicitly 
assumed in NRG studies [18,19] and found also to arise within the LMA [20], this 
remains strictly an assumption (without which essentially no deductions can be made 
about the QCP). 

The spectrum is denoted in this section by D{U;u)) with the {/-dependence 

temporarily explicit (and the rj = 1 + 2ei//7 dependence implicit). As the transition 
is approached, u = \1 — U/Uc{r)\ and the low-energy scale vanishes, as 

= u" (3.1) 

with exponent a. D{U;u}) can now be expressed generally in the scaling form 
ttD{U; u) = m^"''5'q(cj/m") in terms of two exponents a and h (and with a = GFL or 
LM denoting the phase), i.e. as 

7Tu;tD{U;uj) = -^jL.'/iu,) (3.2) 

with the exponent a eliminated and the w-dependence encoded solely in u = w/w*. 
Equation (3.2) simply embodies the universal scaling behaviour of the single-particle 
spectrum close to the transition. Since the QCP must be 'scale-free' (independent 
of w*), it follows from equation (3.2) that the leading x ^ oo behaviour of ^'^(a;) is 
independent of a and given by 

^a{Lo/u.) "/"Ar°° C(r, r;)(|w|M)-'' (3.3) 

with C(r, rj) a constant; such that precisely at the QCP (w* = identically) 

nD{U,;w) = C{r,r])\u>\-K (3.4) 

The exponent b thus governs the w-dependence of the QCP spectrum, and is 
independent of the phase a. In physical terms equation (3.3) embodies the fact that 
the large w asymptotic 'tail' behaviour of the scaling spectra in either phase coincides 
with that of the QCP. More generally of course, the (D-dependence of the scaling 
spectra ^'a(tj) will differ for each phase. 

For the p-h symmetric PAIM (r/ = 0) it is readily seen that the exponent b = r, 
using equation (3.2) and the fact [14] that the leading low-w behaviour of D{U; co) in 
the GFL phase is entirely unrenormalised from the non-interacting limit and given by 

7rD(f/;a;) ''^i^°cos2(|r)|a;|-'' : r] = (3.5) 

We add that the resultant QCP behaviour TrD{Uc;u)) = C(r; 0)|a;|~'' is precisely as 
found in an LMA study of the symmetric PAIM [20], and that the exponent 6 = r is 
likewise consistent with NRG calculations for the GFL phase [16]. 

The behaviour equation (3.5) is exact [14-16], arising because S(a;) vanishes more 
rapidly than the hybridization (oc [^1'"), which itself reflects adiabatic continuity to the 
non-interacting limit. That is manifest also in the limiting low-frequency 'quasiparticle 
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form' for the impurity Green function in the GFL phase. This follows from a low-w 
expansion of the single self-energy 

~ SR(0) - (^1 - 1^ a; (3.6) 

with Z = [1 — {d'£^{u))/dLv)u;=o]~^ the quasiparticle weight which likewise vanishes as 
U Uc{r)—; and with quasiparticle damping embodied in E'(a;) neglected as usual 
[2], as is readily justified perturbativcly. For the symmetric PAIM, where e; = —U/2 
and the impurity charge n = 1, I]^(0) is given trivially by the Hartree contribution 
^Un = \U, whence the renormalised level e* = ei + S^(0) vanishes by symmetry. 
From equations (2.3) and (3.6) the resultant quasiparticle behaviour is then given in 
scaling form by 

7ra;:D(C/;w) ~ 7— , Jj^,-, ^n — ■ V = (3-7) 

* ^ ^ {kiuj + sgn{tj)(3{r)\u)\'-)^ + lu}]^"- ' ^ ' 

from which equation (3.5) is recovered as |d)| — > 0; and where the Kondo scale and 
quasiparticle weight Z are related by 

= [fciZ]T^ (3.8) 

(which relation we add is indeed found within the LMA [15]). 

For the general asymmetric PAIM {r\ ^ 0) we know of no simple argument to 
infer the exponent 6, but NRG calculations arc again consistent with b = r [18] (as 
one might expect from continuity to the symmetric limit). For the asymmetric model 
the renormalised level e* will in general be non-vanishing, and the quasiparticle form 
appropriate to the GFL phase is then given by 

la)!'" 

(fciw - {e*/uji) + iigii(LU)p(r)\Lu\'^y + 

Granted 6 = r it follows from this that the renormalised level e?" must also vanish as 
U Uc{r)— and must do so at least as rapidly as ul, i.e. 

|en~w"' (3.10) 

with exponent cii > ar. The GFL scaling spectrum ^'GFL(a)) = t^^ID{U:ll!) is of 
course obtained by considering in general finite w = w/w*, in the formal limit — > 0. 
From equation (3.9) there are then two a priori possibilities for its \ow-Q behaviour, 
(a) Either ai > ar and thus f* /col 0; whence from equation (3.9) the leading low-d; 
behaviour of the asymmetric GFL scaling spectrum would coincide with that of the 
symmetric limit equation (3.7), namely loID{U;ijj) oc \oj\~^ . Or (b) ai = ar and the 
renormalised level vanishes precisely as col, 

u;. = [\e*\/k]i (3.11) 

with k a constant. In this case the leading low-w behaviour of the scaling spectrum 
is ^'gfl(w) C)c Id'l'', and the quasiparticle form equation (3.9) implies an asymmetric 
double-peak spectrum straddling the Fermi level. It is precisely this characteristic 
behaviour that is observed in NRG calculations [18] (see also §6.1 and figure 11 below), 
from which we conclude that and the renormalised level should indeed be related 
by equation (3.11). NRG calculations for | < r < 1 [18] are moreover consistent with 
an exponent of a = i for the vanishing of the low-energy scale uj^ (oc T* in the notation 
of [18]); implying in turn that the exponent ai = 1. The behaviour just inferred will 
be considered within the LMA in §5,6. We would however add that it gives by no 
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means the whole story, for example the quasiparticle form is by definition confined 
to asymptotically low frequencies Cj and essentially nothing can be deduced from it 
regarding the general form of \1/gfl('^). 

The latter remarks have focused on the GFL phase. But the renormalised level e* 
will in general be non-zero in the LM phase as well, with the leading low-uj behaviour 
of the spectrum thus given from equation (2.3) by nD{U;Lu) oc |w|'"/|ej*p (assuming 
only that T,{ui) vanishes no less rapidly than the hybridization as |w| 0). For this 
behaviour to be part of the scaling spectrum in the LM regime, and because the 
exponent b (= r) of equation (3.2) is independent of the phase, it follows that the 
low-energy scale is again related to the renormalised level by equation (3.11); and 
hence that |e*| must vanish as the transition is approached from either phase. This 
too will be considered within the LMA in §5.1. 

4. Local moment approach 

In this section we develop a local moment approach to the PAIM for arbitrary 
asymmetry 77. For the particular case of p-h symmetry (rj = 0), it reduces to the 
approach given hitherto in [15]. Further details can be found e.g. in [15,27]. 

4.1. Basics and symmetry restoration 

There are three essential elements to the LMA. First that local moments ('/i'), regarded 
as the primary effect of electron interactions [1], are introduced explicitly from the 
outset. The starting point is thus simple static mean-field (MF) i.e. mirestrieted 
Hartree-Fock (HF), containing in general two, degenerate broken symmetry MF states. 
Notwithstanding the severe limitations of MF by itself (see e.g. [15,25,27]), it can 
nonetheless be used as the starting point for a successful many-body theory [25- 
30,15,20]. To that end the LMA employs the two-self-energy description mentioned 
in §2, that follows naturally from the underlying two MF saddle points. The local 
propagator is thus expressed in the form equations (2.9,10), with the two self- 
energies Sct(w) = S(j[{^o-}] (f =T, i) built diagrammatically from the MF propagators 
(themselves denoted by Qr^{ui)). The self-energies are conveniently separated as 

f:,{uj) = t1 + j:^{uj) (4.1) 

into a purely static contribution (E^, figure 1 below) that alone would be retained 
at pure MF level, together with the crucial dynamical piece T,„{ui) that contains in 
particular the correlated electron spin-flip physics intrinsic to both the GFL and LM 
phases. 

For the GFL phase the third key notion behind the LMA, as now discussed 
for the PAIM, is that of symmetry restoration [26-30]: self-consistent restoration of 
the symmetry broken at pure MF level, and hence recovery in particular of the low- 
frequency quasiparticle form (§3) for G{u)) that is symptomatic of adiabatic continuity 
to the non- interacting limit. This is given explicitly by 

(4.2) 

(from which equation (3.9) above follows directly). Equation (4.2) reflects the fact 

that as \uj\ the single self-energy S^(cj) decays to zero more rapidly than the 
hybridization (oc |a;|''), so that neglect of quasiparticle damping is justified and the 
leading low-w behaviour of G{oj) amounts in general to a simple renormalisation of 



— -e* + sgn{uj)l3{r)\uj\'' + isgn(w)|a;|' 
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the non- interacting limit (§3). The question then is under what conditions on the 

two-self-energies {'^^■{uj)} does such behaviour arise? In parallel with the analysis of 
[27] this is readily answered via a simple low-o; expansion of the So-('^)'s, namely 



S^(a;)~S^(0)- 



(J (4.3) 



for the real parts, with = [1 — {dT,^{u)) /du})u,=o\ ^ thus defined; together with 

f:l{uj)^El{u;)r.aM"' (4.4) 
for the imaginary part, with m > r. This form is in fact guaranteed from the diagrams 
(§4.3) for T,„{u) = T,[{Qcr}] with m = 2 + 3r (which concurs with simple perturbation 
theory in U). Equations (4.3) & (4.4) are now used in equations (2.11) & (2.3) to 

determine the resultant \ow-lu behaviour of E(cj) and hence G{lu). From this one finds 
simply that the requisite condition for the quasiparticle form equation (4.2) to arise is 

tf{uj = 0) = Ef (w = 0). (4.5) 

If equation (4.5) is not satisfied, then the resultant leading low-w behaviour of S^(a;) 
is oc \u)\^ which, vanishing as rapidly as the hybridization, vitiates the quasiparticle 
form equation (4.2). 

Equation (4.5) is of course the familiar symmetry restoration condition for 
the GFL phase likewise found for the r — metallic AIM [26,27], and for the 
symmetric PAIM [15] (where it reduces to S^(0) = independently of spin cr, due 
to p-h symmetry). As in previous work [25-30,15,20], it amounts in practice (see 
also §4.2) to a self-consistency condition for the local moment |/i| entering the MP 
propagators {^cr(w)} from which the two self-energies {Ect(i^)} are diagrammatically 
constructed. This generates directly [25,27] a low-energy scale, seen as a sharp 
resonance (at u) = Um.) in the spectrum of transverse spin excitations Im 11"' (to) 
(§4.2); which is non-zero throughout the GFL phase and sets the timescale t ~ h/oj,rn 
for symmetry restoration and hence local 'Kondo singlet' formation. This is the Kondo 
scale, in terms of which single-particle dynamics are found to scale in §6; and with 
Um = = ivk in the notation employed above (low-energy scales are of course in 
general defined only up to an arbitrary constant). 

Imposition of symmetry restoration as a self-consistent constraint thus again 
underlies the LMA for the GFL phase. If equation (4.5) is satisfied then (i) the 
requisite low-energy quasiparticle form is recovered, as above; (ii) all self-energies 
coincide at the Fermi level, namely S^(0) = E^(0) for either a (whence the 
renormalised level equation (2.14) is given equivalently by ef = ei 4- S^(0)); and 
(iii) the quasiparticle weight Z is related to the {Zc} by Z^^ = ^{Z^^ + Z^^). It is 
moreover the limits of solution to the symmetry restoration condition that delineate 
the phase boundaries between the GFL and LM phases (§5.2), as for the symmetric 
PAIM considered in [15] (for U > Uc{r) in the doubly degenerate LM phase itself, 
discussed explicitly in §4.3, symmetry is not of course restored). 

While solution of the symmetry restoration condition equation (4.5) ensures 
correct recovery of the quasiparticle form for the GFL phase, it does not in general lead 
to the Friedel sum rule (equation (2.15)) being satisfied. We have recently developed a 
strategy for metallic AIMs [27] under which both symmetry restoration and the Friedel 
sum rule may be simultaneously satisfied within a two-self-energy framework. That 
approach may be modified straightforwardly to encompass the PAIM, as now sketched 
(and regardless of specific details of the T,a{^^) = ^allQa}] which are considered 
separately in §4.2). 
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Figure 1. Static bubble diagram contribution to So-(t^)- The wavy line denotes 
U, and the solid line the broken symmetry MF propagator. 



4.1.1. Algorithm The two-self-energies are functionals of the MF propagators 
{Gaiuj)} [27], given by 

g^{uj) = [10+ -ei + ax^ M^'T' (4.6) 
where x = |C/|/i| and e; differ in general from their MF values Xq = ^U\fj,o\ and 
eio = 6; + |C/no (with |^ol and no the local moment and charge obtained self- 
consistently at pure MF level). The {Scr(w)} thus naturally depend on both Cj and x: 
Eo-(w) = Scr(w; Ci, x); and, via equations (2.9) & (2.10), so too does G{uj) — G{ui; ej, x) 
and hence Uimp = riimp{ei, x). With such dependence made temporarily explicit, the 
symmetry restoration condition and Priedel sum rule (§2.1) for the GFL phase read 

t,f{0;eux) = t.f{0;eux) (4.7) 

and 

: eliei,x) > 



rii, 



2 : e*{eux)<0 ' ^^'^^ 

These equations are together sufficient to determine both x and ei. In practice, for 
any given U, it is natural to solve equations (4.7) & (4.8) by fixing Ci and determining 
Si, in the following way: (i) For given a solve the symmetry restoration condition 
equation (4.7) for x = ■^U\n\ (and hence (ii) The self-energies {Scr(^^; ei, .x)} arc 
now completely determined, and the bare Ci entering G{u)) (equations (2.9) &(2.10)) is 
then varied until the Friedel sum rule equation (4.8) is satisfied. All quantities are now 
known without looping and the procedure is a rapid numerical exercise. Alternatively 
and equivalently, if one wishes to specify a bare rj = 1 + 2gi/[/ (or ei) from the outset, 
the procedure may be repeated by varying ei until equation (4.8) is satisfied. 

Further, following [27] it is straightforward to show that Ci —a corresponds 
to the particle-hole transformation mentioned in §2, under which Ci — » — [ej -I- U] (i.e. 
77 —7]), nimp — > 2 — nimp and D{uj) D{—uj). As for the metallic AIM [27] we can 
thus focus solely on e; > 0. The case of e; = is of course the particle-hole symmetric 
model [15,25] where the Luttinger integral, equation (2.12), vanishes by symmetry, 
and riimp = 1 for all U; in this case step (ii) above then drops out of the problem and 
the strategy is as presented hitherto in [15]. 

4-2. Self- energies 

We now specify the particular approximation employed for the dynamical contribution 
Sct(w) to the two-self-energies. Equation (4.1) for ^^(a;) is given explicitly by 

t,{uj) = lU{n-a\fi\)+^a{co) (4.9) 
where the first term is the static bubble diagram (figure 1), S° = 
U J_^du) D^^{u;ei,x) with \p\ = |^(ei,a;)| and n = n{ei,x) given by 

l/ij = ^cr / do; D°{u>;ei,x) (4.10a) 

„ J —00 
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Figure 2. Principal contribution to the dynamical £0^(0;). Same notation as 
in figure 1; the transverse spin polarization propagator (see text) is shown as 
hatched. 



(4.10b) 



in terms of the MF spectral density D'^{llj) = —^sgn{Lu)lmQcr{Lo). The class of 
diagrams retained in practice for T,„{co) = I^criiGa}] remains unchanged from previous 
work, and captures the key spin-flip dynamics inherent to both the GFL and LM 
phases. For a full discussion of the diagrams and their physical content the reader is 
referred to [15,25,27]. Shown in figure 2, they translate to 

S,(a;) = r ^ g_„{uj - a;i)n-<^<^(a;i). (4.11) 

Here Il°'~"'{uj) is the transverse spin polarization propagator given by an RPA-like 
particle-hole ladder sum in the transverse spin channel, 

n'"-'"(cj) = °n'^-'^(t^)[i - [/°n'"-'"(w)]-i (4.12) 

where ^Il'^~'^{ui) is the bare particle-hole bubble, itself expressed in terms of the broken 
symmetry MF propagators {Qa} and satisfying the Hilbert transform 

dwi Im°n+~(a;i)sgn(a;i) 



°n+-(a;) 



(4.13) 



We add futher that since [15,27] Tl+-{uj) = (and likewise for the °n'^-'^(a;)), 

only 11"' {ui) need be considered henceforth. 

The functional dependence of the {'Scri'^)} on the MF propagators, for both the 
GFL and LM phases, remains as detailed in [15] for the symmetric model. For the 
GFL phase in particular, 11"' (w) obeys the same Hilbert transform equation (4.13) 
as °n"'"~ and equation (4.11) may thus be written as 



where 



r ^ Imn+-(aa;i) [e{i^^)gz^{uj, + to) + 9i-LUi)gl^iuj, + uj)] (4.14) 

J —CXI 



f 

J — c 



l?O(a;i)0(±u;i) 



a; — Wi ± iO"'" 

denote the one-sided Hilbert transforms such that Qa (w) 



(4.15) 



g+{u^) + g-{uj). 
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Stability The Hilbert transform equation (4.13) for 11+ (w) implies the stabihty 
condition TrRell"' (w = 0) = dwilmll^ (a;i)/|wi| > 0. From equation (4.12), 
together with Im'^II"' {u = 0) = 0, it foUows that this positivity is guaranteed only if 
< U'Re^II"' {co = 0) < 1. An explicit expression for Re°n+~(0) is however readily 
obtained [15,27], 

C/Re°n+-(w = 0) = M^^^ (4.16) 

ImI 

with \jl{ei,x)\ given by equation (4.10a). The pure MF local moment, denoted by 
lHol = |/io(ei)|, is itself given as the self-consistent solution of = \p,{ei,x = ^U\iJ,\)\. 
Equation (4.16) thus shows that the requisite stability is preserved for [/i] > |/j,ol, and 
that the boundary to stability (when t/Rell"' (0; ei,x) = 1) is precisely the pure MF 
local moment |/uo|. 

When = |/io(ei)| it follows from equation (4.16) that Imll^ (w) (equation 
(4.12)) contains a zero-frequency pole. This is of course the essential signature of the 
doubly degenerate LM phase, considered explicitly in the following section. In the 
GFL phase by contrast, = |/io| is determined by symmetry restoration (step 

(i) of §4.1.1); with S\ii\ > correctly found in practice, as required for stability. In 
consequence Imll"' (ui) contains a sharp resonance occurring at a non-zero frequency 
(jj = aJm(= ^*), the low-energy Kondo scale for the GFL phase. This is all entirely in 
keeping with the general picture presented hitherto for the symmetric PAIM, see e.g. 
figure 4 of [15]. As U is increased for any fixed asymmetry r] in the GFL phase, Um 
progressively diminishes ((5|/i| 0) and vanishes at a critical U = Uc{r, rj), where the 
resonance becomes an isolated pole at w = ((5|/x| = 0). This signals the transition to 
the LM phase, to which we now turn. 



4.3. LM phase 

Throughout the LM phase = |/xo(ei)| and hence C/Ren+~(0) = 1. In consequence 
11"' (u) contains an a; = spin-flip pole (of weight Q), together with a continuum 
contribution ^11"' (w) [15], specifically 

n+-(c^) = + 'n+-(a;) (4.17) 

with Q = [U^{dRc"U+- {uj)/du;)^=or^ > 0. lmSn+-(cc>) is given by equation (4.12) 
with 11"' ^11"* , and the real/imaginary parts of ^11^ (w) are again related by the 
Hilbert transform equation (4.13). The basic form for Eo-('^) in the LM phase then 
follows using equation (4.17) in equation (4.11), namely 

E,(a;) = QU^gZ^iio) +^^,(0;). (4.18) 

where ^'E„{u}) is given explicitly by equation (4.14) with 11"' ^H"*"". 

For any given U and 77 (= l+2ei/U) in the LM phase U > Uc{r, rj), Ci is determined 
by requiring that the change in number of electrons due to addition of the impurity 
{nimp, equation (2.6)) is integral, as for the GFL phase discussed in §2.1. Here however 
we require that nimp = 1, i-e. 



2 r 

= -Im / du G{uj) 



1 - 



duj 



= 1 (4.19) 



in common with the local moment regime of the atomic limit (ei < and ei + U > 0) 
to which we expect the LM phase generally to be connected perturbatively in V. 
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Satisfying equation (4.19) within the approach presented thus determines Cj in a 

manner directly analogous to use of the Luttinger theorem to determine e; for the 
GFL phase as described in §4.1.1 (and in fact nimp = 1 is found in practice to be the 
only possible integer solution for nimp in the LM phase). 

Symmetry is not of course restored in the doubly degenerate LM phase, i.e. 
E^(w = 0) ^ t^{Lo = 0). The difference E^(0) - Sj^(O) is however found to be 

sign definite for all U > Uc, such that as U is decreased towards Uc, I'^^iO) — ^^{0)\ 
progressively diminishes and vanishes as the transition is approached, U Uc+- In 
this way the phase boundary between LM and GFL phases may be established, now 
coming from the LM side of the transition U > Uc- And the resultant Uc = Uc{r,r]) 
is correctly found to coincide precisely with that obtained coming from the GFL side, 
U < Uc, as the limit of solution of the symmetry restoration condition equation (4.5) 
for the GFL phase (as indeed is to be expected since as J7 — » Uc+, the resultant 
vanishing of |E|'(0) — E|'(0)| corresponds to the 'onset' of symmetry restoration). 

5. Scales, exponents and phase boundaries 

We turn now to results arising from the LMA specified above. Single-particle dynamics 
will be analysed explicitly in §6fF. Here the behaviour of the relevant low-energy scales 
in the vicinity of the quantum critical point is first considered, followed in §5.2 by 
resultant phase diagrams in the {U,r]) (or equivalently {U, ei)) plane. 



Our aim here is simply to establish the critical behaviour of the relevant low- 
energy scales, specifically cj* (in terms of which single-particle dynamics exhibit 
universal scaling (§s 3,6,7)); and, relatedly, the behaviour of the renormalised level 
e* = €i + E^(0) (equation (2.14)). Unless stated to the contrary, all results given are 
representative of the full range < r < 1 considered in the paper. 

We begin by considering the GFL phase, U < Uc- The critical behaviour of 
the low-energy Kondo scale Wm = is illustrated in figure 3 for r = 0.3 and an 
asymmetry rj = 1 — 2\ei\/U = 0.2; the critical Uc{r, rj) 2± 12.2 being marked as a solid 
square. Specifically, {(jJm/Uy vs U (open circles) is shown, and seen to vanish linearly 
as U Uc—- The critical exponent a of §3 (w* ~ u" with u = \1 — U/Uc\) is thus 
a = 1/r, which agrees with the NRG calculations of [18] (the result is also derivable 
analytically for r <C 1, see §7.1). We find also that ui^ = Wm and the quasiparticle 
weight Z are indeed related by equation (3.8), namely that 



As argued in §3, the renormalised level ej* should likewise vanish as the transition 
is approached, from either phase. When considering the LM phase in particular it will 
in general prove helpful to define spin-dependent renormalised levels via 



by equation (2.11), from which the {e*^} and the renormalised level = Ci -I- E^(0) 
are related by 



5.1. Scales 





(5.21) 
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Figure 3. Critical behaviour of the low-energy scale ui» in the vicinity of the 
QCP; illustrated for r = 0.3 and asymmetry = 0.2 (critical Uc — 12.2, solid 
square), ui^ vs U is shown: for U < Uc (GFL phase, uit, = ujni), [LOm/UY U is 
plotted (open circles); for U > Uc (LM phase, = uij^), |S^(0) - t,f(0)\ oc ojI 

is plotted (open squares). In either ease, ui* ~ |i7 — i7c|~ as the transition is 
approached. See text for further details. 



In the GFL phase (where symmetry is restored, S|'(0) = E|'(0)) the renormahsed 
levels all coincide, ej* = ef^. For the LM phase by contrast, e*,e*| and e^^ are in 
general distinct. 

Figure 4, again for r = 0.3 and rj = 0.2, shows the resultant [/-dependence of e* 
(open circles for U < Uc in the GFL phase, squares for U > Uc). The renormalised 
level is indeed seen to vanish as the transition is approached from cither phase, with 
an exponent ai = 1 (|e*| ^ agreement with the inferences drawn in §3. For 

the GFL phase in particular, and |e*| are thus related by equation (3.11), 

|e*| = kiul (and where the constant k = k(i]) is an odd function of the asymmetry ry, 
such that e* = by symmetry (§3) for the GFL phase of the particlc-holc symmetric 
PAIM, 7] = 0). For the LM phase U > Uc, figure 4 also shows the spin-dependent 
renormalised levels e*^ (triangles); each being seen to be sign definite and to vanish 
linearly as [/ Uc+- 

For the GFL phase, the low-energy (Kondo) scale cj* = arises as a direct 
consequence of symmetry restoration. The reader will however notice that the 
corresponding scale for the LM phase, in terms of which single-particle dynamics 
exhibit scaling and denoted by = wlj has not yet been specified. Here wc 
simply assert, and consider further in §s 6ff, that is related to the spin-dependent 
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Figure 4. Renormalised level e? = ej + S^(0) vs U in the vicinity of the QCP for 
r = 0.3 and r] = 0.2. GFL phase points {U < Uq) are shown as circles, LM phase 
as squares. In the LM phase the t / i-spin renormalised levels e*^ = e; + S^(0) 
are also shown, as up/down triangles. All vanish linearly as the transition is 
approached. 



renormalised levels by 

= WL oc lefi- - el^\-r = |Sf (0) - Sf (0)|^ (5.22) 

(or equivalcntly by cjl oc | e*^. | ; which equivalence follows since each e*^ vanishes 
linearly asU ^ Uc+, figure 4, and low-energy scales are equivalent up to an arbitrary 
constant. The specific proportionality employed in practice will also be specified in 
§6.2). Figure 3 also shows |S|^(0) — Sj^(0)| oc cj^ for U > Uc (open squares); the 
overall figure thus showing the i7-dependence of col in both phases. 

5.2. Phase diagrams 

We now consider representative phase diagrams arising from the LMA, obtained as 
described in §4. As usual the range < r < 1 is considered; and we show in particular 
that the essential characteristics of the phase diagrams divide into two classes, namely 
r < ^ and r > ^. 

Figure 5 shows behaviour typical of r < ^, the phase boundary in the {U, ei)-plane 
for a fixed r = 0.2. As ei is progressively increased, the critical Uc = Uc{r,ei) (solid 
line) first decreases, reaches a minimum at the particle- hole symmetric line Ci = —■^U 
(shown as a dashed line), before rising sharply at Cj ~ —5 and eventually turning 
back on itself. This has two interesting consequences. First that the LM phase is 
inaccessible by increasing U for any fixed ^ — 5 in the present example. Second, for 
fixed ei and due to the above 'turnback', there is at large U a re-entrant transition 
back to the GFL phase (itself seen most clearly in figure 18 below, where a low-r 
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Figure 5. LMA phase boundary in the {U, ei)-plane for a fixed r = 0.2 (solid 
line), representative of the general behaviour for r < i . The dashed line marks the 
paxticle-hole symmetric case: U = —26;. The corresponding MF phase boundary 
is shown as a dotted line. 



phase boundary determined analytically in §7.1 is shown). Both observations are in 
general accord with results of recent NRG calculations [36]; and indeed are seen even 
at the pure MF level of unrestricted Hartree-Fock, shown in figure 5 as a dotted line 
(although the LM phase is strongly exaggerated by pure MF, the many qualitative 
deficiencies of which are discussed e.g. in [15,25,27]). 

A rather more revealing portrayal of the phase boundary for fixed r is shown in 
figure 6 where, for r ~ 0.2, Uc = Uc{r,ri) is displayed as a function of the asymmetry 
T]. Specifically, 77 vs (8/7rr)J7^~^(r, 77) is shown; the constant S/nr simply serving as a 
convenient 'normalisation', it being known from previous LMA work [15] that as r ^ 0, 
U^~^{r,0) ~ 7rr/8 for the symmetric PAIM rj = (which result is exact [15,16]). In 
contrast to Uc{r, Cj) shown in figure 5, which is not single-valued in a (leading to 
re-entrance as above), figure 6 shows Uc{r,r]) to be a single- valued function of the 
asymmetry (which is found for all r); and that LM states arc always accessible upon 
increasing U for any fixed asymmetry |r;| < 1. In §7.1 the 77-dependence of Uc{r,rf} 
arising within the LMA will be discussed further and obtained analytically in the 
small-r limit (which result is also shown in figure 6). 

Figure 6 also shows clearly that particle-hole symmetry is optimal for LM 
formation — 77 = corresponding to the lowest critical Uc{r,r]). For the symmetric 
model it is however well known [10-20] that for r > 1/2 the GFL phase docs not 
survive for any non-zero U. It is for this reason that phase diagrams for r ^ 1/2 
are distinct, and figure 7 illustrates the point. For r = 0.7 the main figure (cf figure 
5) shows the critical Uc{r,ei) vs ei; from which in particular the re-entrant behaviour 
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Figure 6. Phase diagram in tlic (J7, r;)-planc for r = 0.2: t] vs {8/nr)U''~^ is 
sliown (solid line). Tlic corresponding r — + result, obtained analytically in §7.1, 
is also shown (dashed line); see §7.1 for full discussion. 
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Figure 7. Main figure: LMA phase boundary in the (C/, ei)-plaiie for r = 0.7, 
respresentative of the general behaviour for r > ^. The dashed line marks the 
symmetric limit, U = —26;; in this case (77 = 0) only LM states are accessible. 
Inset: phase boundary in the ({/, r))-plane, r] vs {S/'Kr)U^~^ being shown. 




Figure 8. r-dependence of the phase boundaries for fixed asymmetry rj: U~ 
vs r for (from top to bottom) rj = 0.1,0.2,0.3,0.4 and 0.5. See text for details. 
Since the wide-band Umit {D = oo) is employed here, the largest r considered in 
practice is r = 0.95. 



upon increasing U for fixed ej is again apparent. And the inset to figure 7 shows the 
?7-dependence of the phase boundary (in direct parallel to figure 6). The qualitative 
differences between the r < 1/2 and r > 1/2 phase boundaries are self-evident. 

The r-dcpendcncc of the phase boundaries is illustrated in figure 8, which shows 
U~^(r, rj) vs r for a range of asymmetries rj. For rj > the resultant phase boundaries 
have a maximum — in qualitative accord with the NRG results of [13] for [/ = oo 
— which moves to lower r with increasing r). For the symmetric case rj = the 
present approach is known [15] to be deficient as r ^ 1/2, in not capturing the slow 
(logarithmic) divergence of U~^{r — > |,0) seen in NRG calculations [12,13,16]; in turn 
reflecting the neglect of charge fluctuations in the present LMA (which are necessarily 
important in this limit, since Uc{r ^,0) 0). We do not however doubt at least 
the qualitative validity of the results shown for 77 > where Uc{r, rf) remains finite for 
all r considered. 

Finally, note from figure 8 that for small r the critical U~^{r,r]) oc r for fixed 
asymmetry, with a coefficient that is largest in the particle-hole symmetric case 
?7 = 0. We will revisit the matter analytically in §7.1, and show that this behaviour 
is intimately related to the exponential asymptotics of the ubiquitously non- vanishing 
Kondo scale characteristic of the r = metallic impurity model. 




Figure 9. GFL phase spectrum for r = 0.3, asymmetry r/ = 0.2 and interaction 
U = 11 (the critical Uc{r,ri) ~ 12.2); 7r£)(u)) vs w. Insets: spectral behaviour on 
the lowest energy scales. See text for discussion. 



6. Single- particle dynamics 

We turn now to LMA results arising for single-particle dynamics, on all relevant energy 
scales and including in particular the many different facets of their scaling behaviour 
as the underlying quantum phase transition is approached. The GFL and LM phases 
are considered separately in §s 6.1 and 6.2 respectively, and dynamics precisely at the 
QCP in §6.3. 

6.1. GFL phase 

A typical GFL phase spectrum is given in figure 9: TrD(oj) vs oj, shown for r = 0.3, 
an asymmetry r] = 1 + 2ei/U = 0.2 and an interaction strength U = 11 close to 
the critical Uc{r,ri) ~ 12.2. Note first that the change in number of electrons due to 
addition of the impurity is correctly riimp = 0, as required by the Freidel sum rule 
equation (4.8) (where e* = ei + S^(0) > for 77 > 0). By contrast, the local impurity 
charge n = 2 dLvD{uj) is n ~ 0.96, close to the Kondo limit whore n ^ 1; showing 
in turn (since rUmp = 0) that essentially one electron is transferred from the conduction 
band to the impurity upon its addition to the host. 

On high energy scales the expected Hubbard satellites arc seen in the spectrum, 
centred on a; ~ e; and Cj + C/; and correctly broadened two-fold over their counterparts 
at pure MF level for the physical reasons discussed in [15,25]. Spectral evolution with 
U is illustrated in figure 10 where, for the same r and 77, the spectra arc shown upon 
decreasing U further from the phase boundary into the GFL phase. With diminishing 
interaction the Hubbard satellites become progressively less pronounced as expected. 
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Figure 10. GFL phase spectral evolution with U, for r = 0.3 and rj = 0.2 
(the critical Uc ~ 12.2): -kD(lo) vs lo for U = 11,10,9,6 and 5 in obvious 
sequence. Inset: hD{lo) vs \lo\ on a logarithmic scale, for both a; > (solid 
line) and u) < (dashed), w = LOt (= LOm) is indicated by a vertical arrow. The 
function cos^(^r)|a;|~'' is also shown (dotted line); see text for full discussion. 



and ultimately absorbed into the central resonance structure. It is the latter that we 

now begin to consider in more detail, focusing first on the lowest energy behaviour. 

The left-hand inset of figure 9 shows the low-energy behaviour of the spectrum on 
a greatly expanded scale. It consists of two pronounced peaks of unequal intensity, one 
either side of the Fermi level lj = Q, with the more intense peak arising for > 0. This 
is quite general, applying for all r; > (where riimp = 0) and r, for U close to J7c(r, rf). 
For ?7 < this situation is reversed: under 77 —77 as discussed in §2, the spectrum is 
simply reflected, D(uj,ri) = D{—lo, —rj) (and rii^p = 0^2 — nimp = 2). The greater 
peak then lies below, and the lesser above, the Fermi level. For any asymmetry rj 0, 
where e* 7^ in the GFL phase, D{u)) has the ultimate low-w behaviour D{u)) oc |a;|''; 
given trivially (e.g. from equation (3.9)) by 

.Di.) ^ (6.1) 

and symmetric in u. This is illustrated in the right inset to figure 9 (and is of course 
in contrast to the symmetric PAIM [12,14-16,18,20] r? = (and ef = 0), for which 
D{w) oc \u)\~^ as |a;| — » 0, equation (3.5)). 

The qualitative features of the preceding paragraph are in agreement with results 
of recent NRG calculations [18]. We now show that they may be understood using 
the simple low-w quasiparticle form for the impurity spectrum, discussed in §3, 
that is characteristic of the GFL phase. Following that however, we point out the 
severe inadequacies of the quasiparticle form in understanding the GFL phase scaling 
spectrum on all relevant frequency scales, and in particular its essential irrelevance 
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in understanding the behaviour of the spectrum at the quantum critical point (which 
itself is considered in §'s 6.3,7.2). 

The quasiparticle spectrum is given in scaling form by equation (3.9), as a function 
of tD = w/w* with = Wm the low-energy Kondo spin-flip scale characteristic of 
the GFL phase (determined via symmetry restoration as described in §4); and with 
ki = u)l~'^/Z (equation (3.8)) a constant of order unity as the transition is approached, 
where and Z separately vanish (§5). The quasiparticle form is of course restricted 
by construction to <C 1 (which regime is now considered); and for present purposes 
equation (3.9) may be replaced with negligible error (as shown below) by 

nujlD{uj) — (6.2) 

[^-sgn(c.)/3(r)p|12 + p|2. 

where (3{r) = tan(^r) as usual (and with k = \el\/u!l likewise finite as C/ — > Uc{r, rj) — , 
as shown explicitly in §5). Equivalently and conveniently, equation (6.2) may be recast 
as 

'^l^'*!^^"^^ ^ [l-sgn(c.)/3(r)|.;'|1^ + kP'- ^^"^^ 

with (j' = w/lej"!^ (and where e* > has been taken, as appropriate to q > 
0). Equations (6.2) or (6.3) capture the low-energy 'two-peak' spectral structure 
illustrated in the insets to figure 9. Equation (6.1) is of course trivially recovered as 
\uj\ ^0; and from equation (6.3) the resultant spectrum has maxima at a; = ±T*, 
where 

T* = [|er|cos(fr)]^. (6.4) 

For any non-zero asymmetry 77, the low-energy scale T* is proportional to a;* and is 
thus equivalent to it (as follows from equation (3.11), noting that k (= k{r,r])) is odd 
in ?7 such that for the particle- hole symmetric model 77 = 0, = as in §3). We 
have here introduced T* simply to make connection to the NRG calculations of [18] 
for the asymmetric PAIM, where the low-energy scale was defined as the location of 
the spectral maxima. 

In figure 11, LMA results for Tr\e*\D{uj) as the transition is approached arc 
given, plotted vs oj/T* for a range of different r's (and noting that in all cases 
the range \u)\/T* < 10 shown corresponds to \w\/uJt, <C 1, since the constant 
k = \e*\/ujl ~ {T*/uj^Y is found in practice to be small compared to unity). We 
emphasise that the figure shows 'full' LMA results; that the spectra are of course 
scaling spectra, obtained as U ^ Uc{r,r])— and as such universally dependent on 
uj/T* (or equivalently uj) alone; and that the exponent b of §3 (equation (3.2)) is 
indeed b = r. As seen from figure 11 the asymmetrically disposed spectral maxima 
are indeed peaked at iv/T* = ±1, and for the particular case of r = 0.3 the figure 
compares explicitly the LMA results to the approximate quasiparticle form equation 
(6.3). This is not perceptibly distinguishable from the full LMA results, showing that 
equation (6.3) is indeed accurate (the same being true for the other r's shown). 

For sufficiently low oj/T* or ui, the GFL phase scaling spectrum is thus explicable 
in terms of the simple approximate quasiparticle form equations (6.2) or (6.3). An 
obvious question then arises: can this form be extrapolated beyond its strict confines, 
to infer the leading large uj/T* (or a)) form of the scaling spectrum which — as 
shown on general grounds in §3 (equations (3.2-4)) — determines the behaviour of the 
spectrum precisely at the quantum critical point where the low-energy scale oc T* = 




Figure 11. Low-frequency GFL scaling spectra: ■jT\e*\D(uj) vs lu/T* (with the 
low-energy scale T* defined in equation (6.4)), shown for r = 0.1 (dotted line), 
0.2 (short dash), 0.3 (solid), 0.4 (long dash) and 0.5 (point dash). For the case of 
r = 0.3, explicit comparison is also made to the approximate quasiparticle form 
equation (6.3); it is however indistinguishable from the full LMA results (which 
were obtained explicitly for r) = 0.2 but for »? ^ are in fact r^independent over 
the u;/T*-range shown, see equation (6.3)). 



0, and the generic form of which is given by equation (3.4)(7rZ?([/c; w) = C{r;r])\co\~'^)'! 
The reader will notice for example that the leading high frequency behaviour predicted 
from equations (6.2,3) is ttlvIDIu) ~ cos^(^r)|(2)|^'': from which the scale drops out 
of either side, resulting in an apparent QCP spectrum that is indeed of the required 
general form equation (3.4), but with a coefficient C{r,rj) =cos^(-|r). As seen from 
equation (3.5) the resultant QCP spectrum, if correct, would in fact coincide precisely 
with the non-interacting limit of the particle-hole symmetric (ry = 0) PAIM. 

The latter conclusions are not in fact correct, but the above discussion is 
nonetheless informative as now shown. Figure 12 gives the full LMA scaling spectrum 
for r = 0.3 and r] = 0.2: ttujID{lj) for w > (solid line) as a function of a) = w/w*, 
on a logarithmic scale encompassing the entire relevant Q range (and with a) = 1 seen 
as a small 'dip' in the spectrum, a well known and harmless artefact of the RPA- 
like form for IT"' employed here [15,25-27]). The figure also shows the approximate 
quasiparticle form equation (6.2) (dashed line), together with its asymptotic behaviour 
nu)lD{oj) ~cos^(^r)|ci;|~'" mentioned above (upper dotted line). The full LMA scaling 
spectrum is well captured by the quasiparticle form up to uj ~ 10~^, below which is 
indeed seen to exist an appreciable 'intermediate' a)-range over which the behaviour 
TrcjlD{u;) ~cos^(-|r)|£l;|~'' clearly arises. For larger frequencies however, the scaling 
spectrum clearly departs from '~cos^(^r)|d}|~'' behaviour; and for oj ^ 1 is seen to 




Figure 12. Full GFL phase scaling spectrum for r = 0.3 and rj = 0.2: ■kll)^D{ll)) 
vs u) = Lo/ujf {u)t = i^m) plotted for > on a logarithmic scale (solid 
line) and encompassing all relevant ui regimes. The approximate quasiparticle 
form equation (6.2) is shown as a dashed line; the upper dotted line shows the 
rj = non-interacting behaviour cos^(^r)|aj|~''. For |aj| 3> 1 the ultimate high- 
frequency behaviour of the scaling spectrum is also oc but with a different 
coefficient C{r,ri) that is r- and ?7-dependent; the lower dotted line shows the 
latter behaviour extrapolated back to lower lu. Full discussion in text. Inset: 
corresponding results for r = 0.7 and r] = 0.4, same notation as main figure. 



cross over into ultimate large-w behaviour that is again oc \uj\~^ — foUowing generally 
from the fact that the exponent 6 = r as explained in §3 — but with a prefactor 
C(r, ?/) that differs from cos^(^r) and is found in practice to be 77-dependent. 

The above behaviour is not particular to the case r = 0.3 illustrated, being found 
to arise throughout the range < r < 1 considered (the inset to figure 12 shows 
corresponding results for r — 0.7 and 77 = 0.4). The asymmetric GFL phase scaling 
spectrum thus contains two regimes of {Col^^ behaviour: that at lower uj (<C 1) 
corresponds to symmetric (rj = 0) non-interacting limit behaviour, while that arising 
for (D 3> 1 determines the QCP spectrum. This 'two \lu\^^ regime' behaviour is in fact 
known [20] to arise also in the pure symmetric PAIM, 77 = (where the cos^(^r)|a)|~'' 
regime persists all the way down to a; — 0, reflecting the exact low-w behaviour 
equation (3.5) arising in that case). Further consideration of the matter will be given 
in §6.3, and the (r, 77)-dependence of the QCP spectrum arising within the LMA will 
be determined analytically for small r in §7.2. 

Finally, we point out that figure 12 shows the scaling spectrum for w > 0, but 
that it is not of course symmetric in to (or uj /T*) as is obvious e.g. from figure 11. The 
sgn(a;)-dependence of the GFL phase spectrum is illustrated further in the inset to 
figure 10 above (for [r, 77] = [0.3,0.2] and U — 11, quite close to Uc — 12.2). Here the 
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Figure 13. LM phase spectral evolution with U, for r = 0.3 and r] = 0.2 
{Uc ~ 12.2); nD{uj) vs ui for U = 50, 20, 17, 14 (clockwise from top left). 



spectrum nD{u) vs \u)\ (i.e. unsealed) is again shown on a log-scale, for both w > 
(solid line) and u < (dashed). For each sgn(a;) the above 'intermediate' behaviour 
~cos^(^r)|cj|^'" arises (shown as a dotted line). Beyond that however, a crossover 
to ~ C{r,ri)\uj\^^ behaviour occurs; with the C{r,rj) dependent not only on r and rj 
but also for ?7 7^ on sgn(a)) (which behaviour we add is not precluded by the general 
considerations of §3, and which is likewise considered further in the following sections). 
Since the spectrum here is shown in unsealed form the non-universal Hubbard satellites 
are also ultimately visible on the figure, likewise asymmetrically disposed about w = 
because rj ^ 0. Wc also add here that while results from the present LMA are in 
general agreement with the NRG calculations of [18], there is one point on which we 
differ. That we find the coefficient C(r, 77) to depend on 77 (albeit fairly weakly), is 
consistent with a line of (ry-dependent) critical fixed points. Numerical RG calculations 
for the pseudogap Kondo model [13] appear by contrast to be consistent for any r] > 
with a single critical fixed point, and hence in effect no 77-dependence in C{r,ri). 



6.2. LM phase 

We turn now to the LM phase, considering first spectral evolution on 'all scales' as U 
is progressively decreased for fixed asymmetry 77, towards the LM/GFL transition at 
Uc{r,ri). This is illustrated in figure 13, again for r = 0.3 and r] = 0.2 (Uc — 12.2), 
although the behaviour shown is representative of the full r-range considered. TrD{uj) 
vs u! is shown, for U = 50, 20, 17 and 14. For the highest U shown, deep in the LM 
phase, the well developed Hubbard satellites dominate the spectrum entirely, although 
the ultimate low-w spectral behaviour vanishes as oc \w\'^ throughout the LM phase 
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Figure 14. LM phase spectrum nD{uj) vs w for r = 0.3, t] = 0.2 and interaetion 
U = 12.3 elose to the critical Uc — 12.2. Left inset: low-energy behaviour on an 
expanded scale. Right inset: nD(Lu) vs on a logarithmic scale, for both a; > 
(solid line) and u) < (dashed). Pull discussion in text. 



as noted in §3. With decreasing U however, a low-energy spectral structure is seen 
to develop in the vicinity of the Fermi level. It becomes increasingly pronounced as 
U Uc+ and as shown below is the precursor, in the LM phase, of the \uj\~'^ divergent 
behaviour characteristic of the QCP itself. 

For U = 12.3, close to Uc, ttD{uj) vs iv is shown in figure 14; the left inset showing 
the low-energy part of the spectrum. The latter is clearly redolent of its counterpart 
for the GFL phase (figure 9), the low-energy behaviour again consisting of two peaks 
of unequal intensity, the more intense peak arising for w > (and vice versa for 
Tj < 0). In contrast to the GFL phase however, the positions of the peak maxima 
are not symmetrically disposed about the Fermi level; reflecting the fact that in the 
LM phase, where symmetry is not restored, the spin-dependent renormalised levels 
e*^ = Ci + S^(cl; = 0) do not coincide (as shown e.g. in figure 4). In addition, and 
again in contrast to the GFL phase, we find within the LMA that while the lowcst-w 
spectral behaviour is oc as noted above, the coefficient thereof depends on sgn(a;). 
The origins of this are readily seen using equations (2.9,10), which give the leading 
low-w behaviour 

-kD{uj) ~ 2 r5n2 (6-5) 

with the LM phase S<t('^) given explicitly within the LMA considered here by equation 
(4.18). It is the first term of the latter that controls the low-w behaviour of ^^(a^) 
which, using equations (4.6,15), is found to be of form ^^(w) oc \ujy with a coefficient 
dependent on sgn(a;) (as shown explicitly for r? = in [15]); and which, via equation 
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Figure 15. Full LM phase scaling spectrum for r = 0.3 and rj = 0.2, obtained 
as (/ ^ Uc{r,T))+: itljID(lj) vs u) = lo/uj^ (a;* = wl) on a logajrithmic scale, for 
both a) > (solid line) and u) < (dashed). 



(6.5), generates corresponding behaviour in D{lo). 

The low-w peaks in figure 14 (left inset) mark the onset of a crossover from the 
asymptotic oc |a;|'' behaviour, to an \ui\~^ regime which persists essentiahy up to non- 
universal frequency scales on the order ~ 1 (= Aq). This is evident from the right 
inset to figure 14 where the spectrum is plotted on a logarithmic scale: the crossover 
from D{uj) oc \ujY to oc \u>\~'^ behaviour (likewise dependent on sgn(a;), albeit more 
weakly so) is clear, as too are the non-universal Hubbard satellites on the highest 
energy scales. 

The spectra shown in figures 13, 14 arc of course displayed on an 'absolute' scale, 
i.e. shown vs lo in units of the hybridization strength Aq = 1. But the LM phase is also 
characterised (§5) by a vanishing low-energy scale = Wl as the LM/GFL transition 
is approached, so the corresponding LM phase scaling spectrum is readily obtained 
in direct parallel to that for the GFL phase: on decreasing U progressively towards 
?7c-l-, steadily decreases, but for fixed r] the spectra again exhibit universal scaling 
in terms of a) = lo/loi, alone (and with non- universal features such as the Hubbard 
satellites thus 'projected out' as usual). As anticipated throughout we find that it 
is indeed 7r(jj££'(a>) vs u> which exhibits scaling (whence, as for the GFL phase, the 
exponent /; of §3 is again found to be 6 = r as required on general grounds); and with 
Wl given by the form equation (5.22). Specifically we take wl = to be given by 

= If [6*1 - 411" = If (0) - Sf (0)]|^ (6.6) 

where our (free) choice of prefactor, f , is chosen to make ready connection to the 
analytic results for small-r obtained in §7. The resulant LM phase scaling spectrum 
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for r = 0.3,7] = 0.2 is shown in figure 15, showing clearly both the ~ behaviour 
at low-(D and the crossover to the ultimate large-w behaviour oc 

6.3. Quantum critical point 

Which brings us to the quantum critical point itself. From the previous discussion we 
would naturally expect GFL and LM phase spectra arbitrarily close to the QCP to be 
indistinguishable for scales |a;| ^ w*. This is evident in figure 16 which shows spectra 
7rZ)(w) vs u) close to and on either side of the transition \Uc{'i' = 0.3, rj = 0.2) ~ 12.2]; 
for U = 12.0 in the GFL phase and U = 12.4 for the LM phase. On higher energy 
scales the spectra are indeed almost indistinguishable — limited only by the practical 
numerical inconvenience of dealing with the arbitrarily small but strictly non- vanishing 
scales in the immediate vicinity of the transition. 
The inset to figure 16 shows the corresponding a; > spectra on a logarithmic 
scale. The leading low-w behaviour ttD{uj) oc \uj\^ is clear, arising for each phase but 
with different coefficients for the; GFL and LM phases (and each such being r- and 
jy-dependent) . The figure also shows the 'kD{u) oc \co\~'^ behaviour arising for cj ^ u)^ 
in either phase; and, most importantly, that the coefficient of this |a;|~'" behaviour is 
identical for each phase — as argued on general gounds in §3. 

As /7 — » Uc{r, r])± and the low-energy scale — > 0, this power law persists 

to lower and lower frequencies, and alone survives at the QCP where = identically 
and the low-energy physics is scale-free. The LMA enables a direct determination 
of the spectrum precisely at the QCP U = Uc{r,r]) (numerically in general, and 
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analytically for small r as in §7). The resultant QCP spectrum for r = 0.3, = 0.2 is 

also shown in figure; from which the characteristic behaviour ttD{Uc', w) = C(r, r])\Lo\~^ 
as |a;| — > is clearly seen (and which in fact persists essentially up to non- universal 
scales u) ~ 0{1)). 

One obvious but important point should also be noted from the above discussion: 
close to the transition, and in either phase, acts as the natural crossover frequency 
scale to the emergence of common QCP behaviour for ::g> in single-particle 
dynamics, while for lower energies the spectral characteristics arc particular to the 
phase considered. This is of course the dynamical analogue, at T = 0, of the common 
schematic pertaining to static/thermodynamic properties at finite-T (see e.g. [21] or 
figure 1 of [19]), showing the crossover to QCP behaviour upon increasing T above 

6.4- Scaling: a closer look 

In the preceding sections we have focussed primarily on the asymmetric PAIM, rj / 0. 
The purpose of this section is three-fold. First, to compare results for ^ to those 
arising in the particle-hole symmetric limit r/ = [15,16,20]. There are of course 
important differences between the two — notably that for 77 = the LM phase alone 
arises for all r > ^ and any U ^ [10-20] — but, that aside, these are not as radical as 
might perhaps be inferred from the literature and the two may be readily understood 
on a common footing. 

In comparing dynamics for ry 7^ and 77 = 0, an obvious but important point 
should first be noted: that we have considered the scaling spectra in terms of cD = w/w* 
and that the low-energy scale is non-vanishing throughout both phases, including 
the particle-hole symmetric limit rj = [15,20]. In the GFL phase for example one 
might have chosen, as in §6.1, to express scaling in terms of w/T*; with T* cx |e*|~ 
given by equation (6.4) such that the low-energy maxima in the GFL phase spectra lie 
at ui/T* = ±1. This is the procedure adopted in [18] and — in general — is entirely 
equivalent to (w/w,) -scaling because T* cx w*. But the proportionality between the 
two, given (see equation (3.11)) from k{r,ri) = ^ l^)"" '^^ V: such that e* 

vanishes throughout the GFL phase in the symmetric limit, while Wh. = '-^m remains 
non-zero. For this reason, scaling in terms of is naturally required to consider the 
symmetric and asymmetric PAIMs on an equivalent footing. 

In comparing the GFL phase scaling spectra for different 77, it is also convenient to 
consider the so-called modified spectral function [14-16,20] J^{u)) =sec'^{^r)TT\ui\'^ D{uj), 
i.e. 

J^{uj) = sec^{^)r)\u}\^ XTTulD{co) (6.7) 

(likewise universally dependent on Q). For 77 = the resultant spectrum is 'pinned' 

at the Fermi level for all U in the GFL phase [14-16], J^{lu = 0) = 1 (as follows from 
equation (3.5)); while for any ^ 0, oc jw)^'' as \ui\ 0. The scale on which 

this distinction between the r? = and r? ^ scaling spectra arises is T* (oc w*) as 
discussed in §6.1 (figure 11). But in practice the ratio T*/uj^, is found within the 
LMA to be small (r*/w, ~ [k{r,r]]~ with k itself small compared to unity), so the 
qualitative differences between r] = and ^ are apparent only on the lowest-cD scales. 
This is seen clearly in the left inset to figure 17a which, for r = 0.3 and asymmetries 
7] = 0,0.2 and 0.5, shows J^{u)) vs Co for \uj\ <C 1. 
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Figure 17. Scaling spectra in both phases for fixed r (= 0.3) and varying 
asymmetry rj. A (upper): GFL phase scahng spectrum ^{uj) =scc'^ {^r)\uj\'^ X 
TrLj'^D(Lij) vs u) = u)/u)t (t<J* = LOn\) for asymmetries r} = (dotted line), 0.2 
(bold solid) and 0.5 (thin solid line). The main figure shows the central 'Kondo 
resonance'. The right inset shows the spectral 'tails' out to \Cj\ = 150, = 
(dotted line), i) = 0.2 (bold solid line for uj > bold dashed line for tj < 0) 
and ri = 0.5 (thin solid [dashed] line for w > [w < 0]). Left inset: scaling 
spectra for ;/ = 0,0.2 and 0.5 on the lowest scales (c/ figure 11). B (lower): LM 
phase scaling spectrum ixll>^D{ui) vs uj = uj/uit, (wt, = lui^) for rj = (dotted line), 
rj = 0.2 (bold solid line) and rj = 0.5 (thin solid line). Inset: spectra shown out 
to |a;| = 10^ for both tZ) > and < 0, same labelling as in figure 17A. 
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By contrast however, the strong similarities between particle- hole symmetric and 
asymmetric spectra in the GFL phase arc also evident in figure 17a. The main figure 
gives !F{uj) on the |a;| ~ 0{1) scale, showing clearly the low-energy 'Kondo resonance' 
in J^{uj) [14-16] (which becomes precisely the usual Kondo resonance of the metallic 
AIM in the r — limit [14-16]). The effects of non-zero rj, while present, do not lead 
to qualitatively different behaviour either on the ui ~ 0{1) scales or in the spectral 
'tails', shown out to tD ~ C(IO^) in the right inset. These comments hold also for the 
LM phase scaling spectra, t:luI^D{lu) vs Co = uj /uji^ being shown in figure 17b (again for 
r = 0.3 and r] = 0,0.2 and 0.5). In this case 7ra;££'(a;) oc as tD — > for both the 
symmetric and asymmetric cases. The spectra for different r} are thus qualitatively 
similar for all a), albeit that the extent of spectral asymmetry naturally increases with 

•n- 

In previous sections we have emphasised in particular the crossover from the 
oc \ujY behaviour of 'kuj^^D{u}) arising at low-a) to the ultimate largc-cj behaviour 
oc \(b\~^ that in turn controls the QCP spectrum. Our second point is simply to note 
that for the GFL phase, the asymptotic approach to the ultimate \Cb\~'^ form — where 
the corresponding J-^{u}) oc x loI^D{u}) thus plateaus to a constant — is quite subtle. 
In figure 17a for example, where the spectral tails of CF{uj) are shown out to \uj\ = 150, 
has not yet reached its plateau value but is still decreasing, albeit slowly; indeed 
from figure 12 it is clear the plateau docs not occur until cD ^ 0(10'^) for the case 
of r = 0.3 shown. The approach to the ultimate large-w behaviour of the scaling 
spectrum is thus slow, and in fact becomes progessively more so with diminishing 
r ^ 0. The nature of this behaviour (which is relevant to the GFL phase alone) 
is known from our recent work on the symmetric PAIM [20] and is modified only 
quantitatively for r] ^ 0: as r ^ 0, the ultimate large- |a)| behaviour arises only for 
\u}\ ^ e~ . For 1 <C <C e~ by contrast the decay of T{uj) towards its plateau value 
is logarithmically slow. This 'intermediate' large-tD regime progressively dominates 
the GFL scaling spectrum out to ever larger-tD scales as r — » 0, and for r = recovers 
precisely the known [26] logarithmic large-J; tails characteristic of the metallic AIM 
(for which -D(w) oc l/ln'^{\oj\) for \uj\ 1). While it is the small-r behaviour to which 
we turn in the final section, the reader is referred to [20] for further details on this 
point. 

As noted above, the ultimate large-w behaviour u)ID{lo) oc \oj\~'^ is not in fact 
reached in the a)-range shown in figure 17a for the GFL phase. In consequence none 
of the behaviour illustrated in the figure survives at the QCP itself. Our final remark 
here is simply to emphasise this point, evident also in figure 16 (inset): the Kondo 
resonance in !F{u)), illustrated clearly figure 17a and characteristic of the GFL phase, 
'collapses' at the QCP where its weight vanishes; as indeed remarked in [19], and 
which the above figures substantiate. 

7. Small r asymptotics 

Our aim now is to obtain analytical results for the critical behaviour of the low- 
energy scales u;* , the; phase; boundaries and single-particle scaling spectra, in the strong 
coupling regime where Uc{r, 77) » 1. This means in effect that we consider here small 
r <C 1, since oc r as r — > (see e.g. figure 8). 

In strong coupling, where the impurity charge n ^ 1, the low-energy behaviour 
of the pseudogap AIM maps under a Schrieffer- Wolff transformation onto the 
corresponding Kondo model [2], with host density of states /)host(t^) = Po\'^\^0{D—\u}\). 
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The exchange coupHng matrix element J = 2F^[l/|ei| + 1/{U — |ei|)], i.e. J = 
8V'^ /[U{1 — ry^)] expressed in terms of U and 77 = 1 — 2|ei|/C/; and the corresponding 
potential scattering strength is K = rjJ. Since the PAIM hybridization Ai(a;) = 
7rl^^Phost(i^) (with Ai(a;) given by equation (2.4)), it follows that 

1 

ttU 1 ~r;2 

from which strong coupling PAIM results can be transcribed directly to those 
applicable to the Kondo model. 

As a straightforward generalisation of [15,20] we first focus on {ECT(a;)} in the 
GFL phase (equation (4.14)), where in strong coupling Scr(w) takes the form 

^Au)^U^gZ{uj + <jLU^). (7.2) 

This follows from the fact that the spectral weight of Imll^ (ui) in strong coupling 
is confined asymptotically to w > with /Q°°(dw/7r)Imn+~(ci;) = 1 (which physically 
reflects saturation of the local moment, 1); and that its resonance is centred 

on the low-energy spin-flip scale on which scales the {G^{ijJ)} are slowly varying. 
From equation (4.14), equation (7.2) thus results. Equation (7.2) with uim = also 
holds for the LM phase, with uj^^ = reflecting as ever the local spin-degeneracy 
inherent to the LM state; as follows from equation (4.18) noting that [15,20] the 
LM poleweight Q — > 1 in strong coupling (and that the ^So-(w) contribution is both 
negligible in intensity and subdominant in frequency [15]). 

Reg^{u!) is given by the one-sided Hilbert transform equation (4.15), and its 
leading low-u behaviour in strong coupling (and hence for r — > 0) is readily obtained 
using the methods of [15,20,27]. The result is 



TTr 



1_ H 

A'' 



O{\oj\) (7.3) 



where A —imn{D,U) arises as a high-energy cutoff [15,20,27] and = 9a{e[) = 
[1 -I- cre|]~2, with e[ = ^e-jU and Ci determined as specified in §s 4.1.1 and 4.3. 
Equations (7.2) and (7.3) for T,^{uj = 0) wiU now be employed to determine the 
critical behaviour of = Wm, and consequent phase boundaries. We also note that 
the introduction of A =min(£', U) means the following analysis encompasses not only 
D ^ U {as exemplified by the wide band limit D = 00 employed explicitly in the full 
calculations of previous sections), but in addition the case U ^ D. As will be seen, the 
essential differences between the two naturally reside only in the precise dependence of 
the low-energy scales and phase boundaries [15,16] on bare model parameters; while for 
obvious physical reasons the scaling spectra for both phases are entirely independent 
of whether D ^ U ov vice versa. 



7.1. Statics: scales and phase boundaries 

We consider first the approach to the transition from the GFL phase, determining 
thereby the critical Uc{r,rf) for small-r and the critical behaviour of the Kondo spin- 
flip scale cjm = uj* as U ^ Uc — . A brief discussion is then given of resultant phase 
boundaries, and the corresponding low-energy scale = characteristic of the LM 
phase. 
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The symmetry restoration condition equation (4.5) is a central element of our 

approach to the GFL phase. Using equation (4.9) it is given in strong coupling (where 
1) by S? (0) - Sf (0) = U; and hence from equations (7.2,3) by 



AX- 



Prom this the critical Uc{r,r]), where Wm — 0, follows directly from 

(X /Uc{r,r])Y^'^{'Kr/4:){J2a9'r)~^j whence vanishes linearly in r as r — > as 
noted in the discussion of figure 8 (§5.1). In the particle-hole symmetric (r/ = 0) limit, 
e( = (for all r and U, see §4) and hence g+ = 1 = g-; the result obtained previously 
in [15] is then recovered, U~^{r,r] = 0) = irr/S as r — > 0, which is exact [15,16]. More 
generally a knowledge of the as r ^ is required. But since U~^(r,r]) cx r as 
above, the e[ and hence go- may be replaced asymptotically by their corresponding 
r = values [27]; these are [/-independent and functions solely of the asymmetry r], 
given explicitly within the present LMA by [27] 

'^ = yB fiv') = h[^ + {l-v')-]. (7.5) 

(This result may also be obtained by detailed consideration of the limit r — > itself, 

analysis of the renormalised level e* = ei-l-S^(O) showing both that the ga given above 
hold as U ^ Uc{r,ri) and that the exponent for the vanishing of |e*| ~ u""^ (equation 
(3.10)) is indeed ai = 1 as in figure 4.) From equations (7.4,5) the critical Uc{r,r]) is 
thus given from 

A'' r^o 7rr (1 - ry^) 



(7.6) 



or equivalently 



mPoJ. ^ (7.7) 

for the critical exchange coupling of the corresponding Kondo model. 

We shall return to this below, but note first that the U Uc{r,r])— critical 
behaviour of the Kondo scale Wm follows in turn from equation (7.4) as 



The exponent for the vanishing of the low-energy GFL scale Wm = is thus 1/r; as 
found in figure 3 where, with A =min(Z), U) = U as appropriate to the wide band limit 
D = oo there employed, {u)*/Uy vs U was shown to vanish linearly in m = {1 — U /Uc) 
as u — > 0+. We add moreover that using equation (7.6) in equation (7.8) and taking 
the limit r — > recovers precisely the Kondo scale hitherto found in [27] for the r = 
metalUc AIM, 



A --"v 8 /(„») ; - "V p„jfm 

(with J = J{U, r]^) given by equation (7.1)). While recovery of an exponentially small 
Kondo scale from an approximate theory is non-trivial, we note that the exponent 
therein differs in general from the exact result [2] for the r = Kondo model by the 
factor of f{rf) G [1, ^]. It is thus as such exact only for the symmetric model where 
f{T] = 0) = 1; albeit that the correction is modest, / being slowly varying in ry and lying 
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Figure 18. Phase boundary in the (?7, ei)-planc for fixed r = 0.05 (solid line), 
obtained from equation (7.6) (with \ = U as appropriate to the wide band case). 
The corresponding result for r = 0.01 is shown for comparison (dashed line). 
Inset: r = 0.05 phase boundary for large negative e;, which is seen to approach 
asymptotically Uc = |ei|+- 



e.g. within 10% of unity for r] < 0.6. Since equation (7.6) (or equivalently equation 
(7.7)) likewise reflects the r — exponent, the 77-dependcncc of the corresponding 
exact result as r — » is obtained by replacing /(?7^) by 1 therein. 

The r ^ phase boundary arising from equation (7.6) is shown in figure 6 of 
§5.2 (with \ = U as appropriate to the wide-band limit), where it is compared to 
full numerical results for r = 0.2 and displayed as rj vs {8/nr)U^~^{r,r]) — which we 
note is equivalently limr^o[Uc{r,r] = 0)/Uc{r,'r])]. Here we add in passing that the 
resultant r ^ phase boundary between non-magnetic (GFL) and magnetic (LM) 
states, is redolent of the phase diagram obtained in Anderson's original paper on the 
r = metallic AIM (figure 4 of [1]). There is of course no transition for r = 0: the 
'phase diagram' of [1] was simply the boundary to local moment formation at MF 
level, the qualitative deficiencies of which are long since known; in particular that a 
spurious transition arises at a critical U^^{r = 0,r]) (with [7^^(0,0) = tt). But the 
similarities between the two arc nonetheless striking and the r ~ MF results for 
ri vs [U^^{0,r] = 0)/C/^^(0, r;)] ([1], figure 4), while wrong for r = itself, appear 
ironically to be qualitatively correct as r — > 0. 

Figure 18 by contrast shows the phase diagram in the (C/, ei)-planc arising from 
equation (7.6) (with A = U), for both r = 0.05 (solid line) and r = 0.01; the former 
showing again the re-entrant behaviour discussed in §5.2 in connection with figure 
5. The increasing stability of the GFL phase with diminishing r is clearly evident, 
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0.2 0.4 0.6 0.8 

r 



Figure 19. Full LMA phase boundaries, siiown as F(r,r)) = (8/7r)[/(»?^ )/(l — 
ri'^)]U^~^ vs r, for different asymmetries rj = 0.1, 0.2, 0.3, 0.4 and 0.5 (from top 
to bottom in sequence). As r —* 0, all collapse to the common form F{r,T]) ~ r 
(itself shown, dashed line). 



and expected physically since the Fermi liquid phase alone persists throtighout the 
entire {U, ei)-plane in the r = limit of the metallic AIM. The inset to the fignre also 
shows the phase boundary (for r = 0.05) out to large negative ei, from which it is seen 
that the critical Uc asymptotically approaches |ei|+ (shown as a dotted line) — again 
physically natural, this being the condition for local moment formation in the atomic 
limit. 

In figure 19 the r-dependence of full LMA phase boundaries is shown for a range 
of asymmetries ri, and the wide-band limit employed in previous sections. Specifically, 
F{r,r]) = (8/7r)[/(?7^)/(l — r]'^)]UJ^~^{r,T]) vs r is shown, which from equation (7.6) 
should behave as F{r,r]) ~ r as r ^ such that phase boundaries for difi'erent r 
then collapse to a common curve. This behaviour is indeed seen in figure 19 and 
we note that, in line with results from a previous NRG study [13], the resultant 
phase boundaries have little asymmetry dependence up to r ~ 0.25. The asymptotic 
behaviour F{r,r]) ~ r is also shown on the figure and seen to be reached in practice 
below r ^ 0.1 or so (which range we add is known from a previous LMA and NRG 
study [16] of the symmetric limit to be enhanced further if one works instead with a 
finite bandwith D, with U/D > 1). 

We now consider briefiy the corresponding U Uc+ critical behaviour of the 
low-energy scale = characteristic of the LM phase, given by equation (6.6) 
(the resultant phase boundary Uc (r, ij) is of course correctly independent of the phase 
from which it is approached) . This is straightforward, since in strong coupling (where 
ImI ^ 1), tf{0) - tf{0) = (0) - T,f{0) - U from equation (4.9); and S^(0) 
is given from equations (7.2,3) with u>m = throughout the LM phase, whence 
Sf (0) - Ef (0) = {Ay/nr) J^a 9<t - U. But from equation (7.4) {AY /-Kr) 9<t = Uc, 

so (0) - Sf (0) Uc(U/Uc - 1) as t/ ^ Uc+, with C/c oc l/r as r ^ from 

equation (7.6). Hence wl oc (r|S?'(0) — E|'(0)|)^/'' given by equation (6.6) indeed has 
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the critical behaviour 

with exponent 1/r, as shown numerically in figure 3 (§5.1). 
7.2. Scaling spectra and the QCP 

We turn finally to the scaling spectra characteristic of the GFL and LM phases (figures 

12, 15), in particular the leading large uj = uj/uj^, behaviour thereof that determines 
the spectrum at the quantum critical point itself (§6.3). 

The impurity Green function is given by equations (2.9,10), which may be 
expressed quite generally as G{uj) = IX^ctI^^ ^ ^{^) ^ ~ (^^(w) — T,c{0))]~^, 
where we use Scr(u;) — S(j(0) = So-(w) — Scr(O) (trivially, from equation (4.9)). Here 
as usual, e*^ = ei + 2^(0) is the renormalised level (equation (5.20)); with the e*^ = ej* 
independent of spin cr for the GFL phase where symmetry is restored (equation (4.5)), 
but ^ e* for the LM phase, as shown in §5.1. As discussed in the preceding sections, 
(i) the scaling spectra for the GFL or LM phases are obtained by considering finite 
w = w/w* (with (J* = Wm or wl as appropriate), in the formal limit 0; and (ii) it is 

u)lG{ijj) (as opposed to G{u!) itself) which exhibits universal scaling in terms of uj. This 
follows directly tolGico) ^ i EM""^^ - <^7''M^) - ^^{^1 + (S.(co) - S,(0))]-i, 
in which a'^~''|ci)| may be neglected entirely for < r < 1 as considered. By contrast, 
from equations (2.4,5) for the hybridization, uj~'^A{uj) is non-vanishing in the scaling 
regime and given by w^^A(w) = — sgn(a;)[/?(r)-|-i]|tD|'" (independently of the bandwith 
D). The scaling spectrum 'jtu)ID{u) = —sgn{u)wllmG{u)) is thus obtained from 

-1 



sgn(o.)[/3(r) +i]|c:.r " ^ " <'{^Au^) " S.(0)) 



(7.11) 



which we add holds generally for all < r < 1 considered; and where e^^/ujl as 
^ has been shown in §5.1 to tend to a non-zero constant for both the GFL and 
LM phases (denoted therein as k{r,ri) for the GFL phase). 

Since ulGioj) scales, it follows from equation (7.11) that uj^^'lYiaiuj) — Scr(O)] 
is likewise universally dependent solely on uj — lu/lo,:. In the strong coupling, r — > 
regime of interest here, this may be shown directly using equations (7.2,3) (together for 
Ej,(w) with Img^iuj) = -sgn{uj)TrD"jLj)e{±uj) with D°{uj) the MF spectral density). 
For the GFL phase one thereby finds 

uri^^iu) - E^(0)] = -a^g, [\co + - 1] (7.12a) 

w-'-S^(w) = 4g, e{-[l + aoj])\u> + a]'' (7.12b) 

(noting that I]J^(0) = 0); where the = gcr{r]) are given explicitly by equation (7.5) 
and = Wni- For the LM phase by contrast, 

coZ^i^^iu;) - S^(0)] = -a—g^ (7.13a) 

LJ-^Ki^) = 45a e(-aw)|a;r (7.13b) 

with Lu^ = t^L. Equations (7.12,13), which we emphasise hold for all uj (and r ^ 0), are 
indeed of the required universal form; and for 77 = 0, where the ga = I, reduce to the 
results of [20] for the symmetric PAIM. Together with equation (7.11) they determine 
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the full (D-dependence of the scaling spectra for the respective phases. The resultant 

range of behaviour is rich, encompassing for example the intermediate logarithmic 
regime in the GFL phase mentioned in §6.4, and may again be obtained following the 
arguments of [20]. 

Here we consider solely the ultimate large-(D behaviour of the scaling spectra for 
either phase, relevant to the QCP itself. This arises for \u!\^ 1, and for r ^ (where 
/3(r) =tan(|r) ^ 0) equations (7.11-13) yield the requisite behaviour 



given explicitly for positive uj (for < 0, g+ and g_ arc simply interchanged). The 
result equation (7.14) applies to both the GFL and LM phases; and cancels out 
of each side of the equation to give the behaviour 7tD{u>) = C{r,r])\Lj\~'^ precisely at 
the QCP (which in that case extends down to uj = (§6.3)). As foimd numerically 
the QCP coefficient C{r,r]) is indeed both r- and Ty-dependent; with C(r,r]) ~ as 
r — > and as such from equation (7.6) proportional to l/[J7c(r, rj)]^ (or equivalently 
[PqJc^)i in which the interacting nature of the QCP is directly manifest. 

8. Summary 

We have developed in this paper a non-perturbative local moment approach to the 
T = pseudogap Anderson impurity model, with a conduction electron density of 
states proportional to \lo\^ and itself a longstanding paradigm [6-20] for understanding 
the topical and important issue of local quantum phase transitions. The transition 
here is between a generalised fermi liquid phase, perturbatively connected to the non- 
interacting limit, in which the impurity spectrum exhibits a Kondo-like resonance 
indicative of underlying local singlet behaviour; and a doubly degenerate local 
moment phase in which the impurity spin remains unquenched. In addition to phase 
boundaries, and associated critical behaviour of relevant low-energy scales, local single- 
particle dynamics in the immediate vicinity of the transition have been considered: in 
the GFL and LM phases separately, and at the quantum critical point itself where the 
Kondo resonance has just collapsed and the problem is scale-free. Notwithstanding 
the range of vanishing low-energy scales as the transition is approached, and their 
importance in understanding local spectra, dynamics close to the QCP in either phase 
arc nonetheless controlled by a single low-energy scale cj, and as such scale universally 
in terms of w/w* . A determination and understanding of these scaling dynamics on 
all w/cj* scales, including their continuous evolution 'from the tails' to pure QCP 
behaviour, has been a primary focus of the work. Approximate though it is the LMA 
is to our knowledge one of the few theoretical approaches currently capable of handling 
these issues; and the description arising from it is both rich and arguably a good deal 
more subtle than appears to have been appreciated hitherto. Much of course remains 
to be understood, not least the question of dynamics at finite-temperature, aspects of 
which will be considered in a forthcoming paper. 
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